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AGNPS,  AGRICULTURAL  NON-POINT-SOURCE 
POLLUTION  MODEL 

A  WATERSHED  ANALYSIS  TOOL 


Robert  A.  Young,  Charles  A.  Onstad, 
David  D.  Bosch,  and  Wayne  P.  Anderson 

INTRODUCTION 

The  importance  of  runoff  from  agricultural 
lands  as  a  nonpoint  source  of  pollution 
has  brought  about  an  effort  in  Minnesota 
to  develop  a  uniform  method  of  analyzing 
the  quality  of  runoff  from  agricultural 
watersheds.   Federal  Law  92-500,  Section 
208,  has  required  that  all  States  evaluate 
upland  erosion  and  determine  its  effect  on 
water  quality.   As  part  of  this  effort, 
objectively  evaluating  potential  pollution 
problems  on  agricultural  watersheds  within 
the  State  of  Minnesota  is  necessary. 

A  few  years  ago  the  Minnesota  Pollution 
Control  Agency  (MPCA)  realized  the  need 
for  a  uniform  method  for  evaluating  agri- 
cultural watersheds  in  the  State.   The 
MPCA,  the  Minnesota  Soil  and  Water  Conser- 
vation Board,  and  the  U.S.  Department  of 
Agriculture's  Soil  Conservation  Service 
and  Agricultural  Research  Service  entered 
into  a  memorandum  of  understanding  to 
develop  a  model  that  could  analyze  both 
sediment  and  nutrient  transport  in  a 
watershed.   The  model  was  to  be  developed 
for  a  large  computer  system  to  analyze 
agricultural  watersheds.   The  objectives 
of  the  model  were  to  obtain  uniform  and 
accurate  estimates  of  runoff  quality  with 
primary  emphasis  on  sediment  and  nutrients, 
to  compare  the  effects  of  various  conserva- 
tion alternatives  on  implementation  as  part 
of  the  management  practices  of  the  water- 
shed, and  to  develop  a  flexible  and  easy- 
to-use  model. 

The  Agricultural  Non-Point-Source  Pollu- 
tion Model  (AGNPS)  presents  the  user  with 
a  means  of  objectively  evaluating  non- 
point-source  pollution  from  agricultural 
watersheds.   The  input  data  for  the  model 
can  be  obtained  from  readily  available 


records  and  visual  reconnaissance  in  about 
1  person  month  for  larger  watersheds  (up 
to  23,000  acres  in  size),  and  about  3 
person  days  for  smaller  watersheds  (up  to 
500  acres  in  size).   Analysis  of  pollutant 
loads  from  feedlots,  investigation  into 
the  effects  of  implementing  various 
conservation  practices  including  impound- 
ment terraces,  and  the  ability  to  output 
water-quality  characteristics  at  interme- 
diate points  throughout  the  watershed 
network  are  all  within  the  model's 
capabilities.   AGNPS  is  a  single-event- 
based  model  intended  to  simulate  sediment 
and  nutrient  transport  primarily  from 
agricultural  watersheds  in  the  State  of 
Minnesota,  although  the  principles  on 
which  the  model  is  based  are  not  limited 
to  Minnesota.   Proceeding  from  the 
headwaters  of  the  watershed  to  the  outlet, 
pollutants  are  routed  in  a  stepwise 
fashion  so  the  flow  at  any  point  may  be 
examined. 

The  model  works  on  a  cell  basis.   These 
cells  are  uniform  square  areas  that  divide 
up  the  watershed.   This  division  makes  it 
possible  to  analyze  any  area  in  the 
watershed.   The  basic  components  of  the 
model  are  hydrology;  erosion;  sediment 
transport;  and  transport  of  nitrogen, 
phosphorus,  and  chemical  oxygen  demand. 
In  the  hydrology  portion  of  the  model, 
calculations  are  made  for  runoff  volume 
and  peak  concentrated  flow.   Total  upland 
erosion,  total  channel  erosion,  and  a 
breakdown  of  these  two  sources  into  five 
particle  size  classes  (clay,  silt,  small 
aggregates,  large  aggregates,  and  sand) 
for  each  of  the  cells  are  calculated  in 
the  erosion  portion.   Sediment  transport 
is  also  calculated  for  each  of  the  cells, 
in  the  five  particle  classes  as  well  as 
the  total.   The  pollutant  transport  por- 
tion is  subdivided  into  one  part  handling 
soluble  pollutants  and  another  part 
handling  sediment-attached  pollutants. 


USER'S  GUIDE  TO  AGNPS 


Introduction  to  User's  Guide 

This  section  explains  how  to  collect  and 
compile  information  needed  to  run  AGNPS. 
It  will  also  serve  the  computer  programmer 
as  a  reference  in  the  establishment  of 
data  files. 

This  section  outlines  and  explains  the 
variables  needed  for  the  model  and 
describes  where  they  may  be  obtained.   It 
also  explains  and  provides  examples  of  two 
methods  of  establishing  input  files.   The 
first  method  is  the  manual  collection 
method.   In  this  particular  method,  the 
user  collects  all  the  variables  needed 
for  running  the  model  without  aid  from  an 
outside  data  source.   The  second  method 
makes  use  of  the  Minnesota  Land  Management 
Information  System  (MLMIS),!  a  geographic 
information  system  that  contains  detailed 
data  for  many  parts  of  the  State.   Much  of 
this  information  fits  the  needs  of  the 
model  and  can  be  used  to  establish  a  part 
of  the  input  file. 

AGNPS  was  written  in  the  FORTRAN  IV 
computer  language  and  developed  on  a 
Hewlett-Packard  1000  computer  system. 
Further  technical  information  on  the  model 
may  be  found  in  the  second  section  of  this 
book,  "Model  Documentation,"  which 
documents  the  methods  and  equations  used 
in  the  model. 

Model  Basics 

AGNPS  is  intended  to  provide  basic 
information  on  water  quality  to  be  used  to 
classify  non-point -source  pollution 
problems  in  agricultural  watersheds.   The 
model  provides  outputs  on  hydrology, 
with  estimates  of  both  volume  and  peak 
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runoff,  and  on  sediment,  with  estimates  of 
upland  erosion,  channel  erosion,  and 
sediment  yield.   Along  with  these,  the 
user  will  receive  estimates  of  the 
pollutants  nitrogen  (N) ,  phosphorus  (P), 
and  chemical  oxygen  demand  (COD),  in  units 
of  concentration  and  mass,  contained  in 
the  runoff  and  the  sediment.   The  outputs 
can  be  examined  for  either  a  single  cell 
or  for  the  entire  watershed.   This 
information  may  then  be  used  to  rate  the 
watershed  objectively  against  other 
watersheds,  to  further  pinpoint  water- 
quality  problems,  and  to  investigate 
possible  solutions  to  these  problems. 

Watersheds  examined  by  the  model  must  be 
segmented  into  square  working  areas  called 
cells.   This  allows  the  model  to  provide 
water-quality  information  for  specific 
locations  within  the  analyzed  watersheds. 
Groupings  of  the  cells,  determined  by 
drainage  pattern,  make  up  numerous  sub- 
watersheds.   These  subwatersheds  may  be 
separately  analyzed  to  provide  further 
information  on  individual  watershed 
sections . 

Data  required  for  model  execution  may  be 
obtained  through  MLMIS,  visual  analysis, 
maps  (both  topographic  and  soils)  and 
various  technical  publications,  tables, 
and  graphs,  either  included  in  this  manual 
or  easily  acquired  elsewhere.   Data  can  be 
classified  into  two  categories:   watershed 
data  and  cell  data.   Watershed  data 
include  information  applying  to  the  entire 
watershed  and  to  the  storm  event  to  be 
simulated.   Cell  data  include  physical 
information  describing  each  of  the  cells 
as  well  as  information  based  on  the  land 
practices  in  the  cell.   A  numbering  system 
is  used  to  label  each  of  the  cells  and  to 
organize  the  collected  data.   Sample  data 
collection  forms  are  illustrated  in 
figures  1-6. 
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Figure  2 

Watershed  data  collection  sheet. 
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Figure    3 

Feedlot  data  collection  sheet. 
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Establishing  a  Watershed  Data  File 


Preliminary  Examination 

A  preliminary  investigation  of  the 
watershed  is  necessary  before  the  input 
data  file  can  be  established.   The  steps 
to  this  preliminary  examination  are: 

Step  1 

Obtain  a  detailed  topographic  map  of  the 
watershed  to  be  analyzed.  We  recommend 
a  U.S.  Geological  Survey  map,  scale 
1:24,000.   Although  this  particular  map  is 
not  required,  accuracy  and  uniformity 
make  it  desirable. 

Step  2 

Using  the  topographic  lines  of  the  map, 
establish  the  drainage  boundaries  of  the 
watershed.   If  MLMIS  is  to  be  used  as  a 
data  source,  the  boundaries  of  the 
watershed  will  be  documented  by  them.   The 
boundaries  on  their  files  must  match  those 
used  in  the  preliminary  examination  if 
their  data  are  to  be  used. 

Step   3 

Divide    the  watershed    into    cells    as    shown 
in  figure   7.      By   quartering   first    the 
sections    and    then   the   quarters,    you   can 
divide   the  watershed   into   40-acre  cells. 
If   necessary,    the   cells   can  again  be 
quartered    down   to    10   acres.      The 
guidelines    for   the   cell   sizes    are: 


Cell 

dimensions 

(ft) 


Watershed 

Cell 

area 

area 

(acres) 

(acres) 

<2,000 

10 

>2,000 

40 

660  x   660 
1,320  x  1,320 


Only  those  cells  with  more  than  50  percent 
of  their  area  within  the  watershed 
boundaries  should  be  included.   If  the 
MLMIS  data  files  are  to  be  used,  their 
cells  must  match  those  drawn  up  by  the 
user  for  compatibility  in  later  analyses. 


Watershed 
Boundaries 


Figure  7 

Sample  division  of  a  watershed 

into  cells. 


On  receiving  the  MLMIS  file,  lay  out  the 
grid  system  used  on  the  map. 

Step  4 

Number  the  cells  consecutively  from  1  to 
the  number  of  cells,  beginning  at  the 
cell  in  the  northwest  corner  of  the  water- 
shed and  proceeding  from  west  to  east 
southward  (fig.  8).   This  numbering  system 
is  used  in  the  model  for  labeling  cells 
and  will  aid  in  quickly  identifying  the 
cells  in  program  output. 

Step  5 

Establish  the  watershed  drainage  pattern 
from  the  cells.   The  cell  drainage  direc- 
tion is  defined  as  the  direction  of  flow 
leaving  the  cell.   This  can  be  one  of 
eight  possible  directions,  directly  out 


I  /  Watershed 
||    Boundaries 


Figure   8 

An  example   of    the   numbering   system 
used    to   identify   cells    and    the 
drainage   patterns    in  the   cells. 

the   sides    of    the   cells    or  out    of    the 
corners    (fig.    9).      This   step   is  most 
easily   accomplished    by  drawing   an  arrow 
out    of    the    cell   in   the   direction  of    flow, 
as   in  figure  8. 

Forming   a  Data  File 

Once   the   preliminary   examination   is 
complete,    the   input    data  file    can  be 
established.      The   data   can   be   obtained 
from  public    records    and   personal    inspec- 
tion and    through  use    of    the   information  in 
this   book.      An  MLMIS   data   file    can  also   be 
used   to   provide   all   or   part    of    the   infor- 
mation.     The   data   file   parameters    and    the 
required    format    for   input    of    this    file 
into   the   computer  model    are   shown   in   table 


Figure   9 

Identification  numbers   given   to   each 

drainage   path  leaving   the   cell. 


1.      An   explanation  of    the   data  file   and 
information  outlining   the   collection  of 
the   data   follows.      The   steps   required 
repeat   and   expand   on   some   of    the   steps 
used    in   the   preliminary   examinations. 

Line   1 

Watershed  identification:   This  is  a  name, 
30  characters  or  fewer,  given  to  each 
watershed.   It  identifies  the  input  and 
output  data  files  for  the  watershed. 

Line  2 

Area  of  each  cell:   Each  watershed  is 
divided  into  square,  uniform  units  called 
cells  (fig.  7).   The  size  of  these  cells 
depends  on  the  degree  of  detail  desired  in 
the  analyses  and  on  the  size  of  the  water- 
shed.  For  watersheds  up  to  2,000  acres  in 
size,  we  recommend  10-acre  cells,  and  for 
larger  watersheds,  40-acre  cells.  But,  if 
a  very  detailed  investigation  is  desired, 
you  can  use  2.5-acre  cells  for  the  smaller 
watersheds  and  10-acre  cells  for  the 
larger  ones.   Likewise,  you  can  use  larger 
than  recommended  cells. 
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Table    1. 

Data  file   parameter  explanation  and    input   format 


Basic  watershed    (no  impoundment    terraces    or   feedlots)* 
Line  Parameter 

1  Watershed   identification   (description,    30-character  maximum) 

2  Area  of   each   cell   ( acres )/nuraber   of    cells/precipitation  (inches)/ 

energy-intensity  value   of    the   storm/description     Format   (F4.1, 
I4,2F6.1,10X,A) 

3-end  Cell   parameters 

of 
file 

Column      Content 

1  cell  number 

2  receiving   cell  number 

3  SCS  curve   number    (table    2) 

4  land   slope    (percent) 

5  slope   shape   factor:      1 — uniform,    2 — convex,    3 — concave 

6  field  slope    length   (feet)    (table   4) 

7  channel    slope    (percent) 

8  channel   sideslope    (percent) 

9  Manning's  roughness   coefficient    for   the  channel    (table   5) 

10  soil   erodibility   factor 

11  cover   and   management   factor   (appendix   5) 

12  support    practice   factor   (Wischmeier   and   Smith   1978) 

13  surface   condition  constant    (table    2) 

14  aspect    (number   showing   direction  of   drainage,   fig.    9) 

15  soil    texture:      0 — water,    1 — sand,    2 — silt,    3 — clay, 

4 — peat 

16  fertilization  level:      0 — none,    1 — low,    2 — medium,    3 — high 

17  fertilizer  availability   factor   (percentage   of    fertilizer 

left   in  top   one-half    inch   of    soil,   table   7) 

18  point   source   indicator   (a   single-digit  number   other   than 

0  denotes    the  number   of    sources    in  the  cell) 

19  gully  source   level   (tons) 

20  chemical   oxygen  demand    (COD)    factor   (table   8) 

21  impoundment    factor    (number   of    ponds    in   the   impoundment 

terrace   system,   maximum  of    13) 

22  channel    indicator   (a   single-digit  number   other   than  0   denotes 

the  number  of   channels    in  the  cell) 
Format        (3I4,F5. 1 ,I2,I4,2F5. 1 ,F5.3,2F4.2,F5.2,F4.2,3I2,I4,I2, 
214,13,12) 


See   footnote  at  end   of   table. 
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Table      1. — Continued 

Data   file   parameter   explanation  and    input    format 

Watersheds   with  one   or  more   point    sources  * 
Line  Nonfeedlot    point    source   parameters 

1  1/flow   (cfs)/TN   (p/m)/TP(p/m)/COD(p/m)      Format    ( 11 , I5,3F7 . 1 ) 


Feedlot  parameters 


1  2/feedlot  area  (acres)/f eedlot  curve  number   Format  (I1,F6.2,I4) 

2  Acres  of  roofed  area/area  2  subareas  area  (acres)/subareas  curve 

numbers  (up  to  six  subareas)   Format  (7F6.2,6I4) 

3  Area  3  subareas  area/subareas  curve  numbers   Format  (6F6.2,6I4) 

4  Buffer  area  slope/surface  condition  constant/flow  length  (feet) 

(up  to  3  sets  of  data  for  a  buffer  area)   Format  (F5. 1  ,F5.2 ,  14, 
F5.1,F5.2,I4,F5.1,F5.2,I4) 

5  Number  of  animal  type/COD  factor/phosphorus  factor/nitrogen 

factor  (up  to  three  different  animal  types)   Format  (I5,3F5.2, 
I5,3F5.2,15,3F5.2) 


Watershed  with  one  or  more  impoundment  terraces  1 


Line     Parameter 


Area  draining  into  the  impoundment  (acres)/diameter  of  impound- 
ment's pipe  outlet  (inches)  Format  (F4.1,I2)  (up  to  13  pairs 
of  data  in  this  format) 


^If    an  impoundment    terrace   or   a   point    source    such   as    an  animal   feedlot, 
spring,    waste   treatment    plant,    etc.,    is    in   the   cell,    more   information 
is   needed — one   line  of    impoundment   data  or  nonfeedlot    point-source 
data,    if    any,    followed   by   5    lines    of    feedlot    data,    if    any. 


1  2 


Number   of    cells 
the   watershed. 


The   number   of   cells    in 


Precipitation:      To  make   uniform   compari- 
sons  among  watersheds,    you  will  normally 
input   a   precipitation  value   in  inches 
appropriate   for   the   particular  watershed 
area  for   a  24-hour,    25-year   storm  frequen- 
cy.     If    further  analysis   is  desired, 
different    storms   can  also  be   simulated   by 
the   model.      See   appendix   1    for   the   rain- 
fall  frequency  values   for  Minnesota. 

Energy-intensity  value:      This  value    is 
the  rainfall   erosion  index  used   in  the 
universal   soil   loss  equation  (USLE).      The 
value   must  be   for   a   single   storm  event 
only,    in  units  of    foot-tons    per   acre-inch. 
The   energy-intensity  values   appear   along 
with   the  rainfall   frequency  curves   in 
appendix   1. 

Description  (optional):      A  brief    statement 
describing   specific   conditions    that   the 
data  represents,    i.e.,   year,    conservation 
practices   installed,   etc. 

Line  3   to   end   of   file 

Cell  Parameters 

1.  Cell  numbers:      Each   cell   in  the 
watershed   is   identified   by  a  number. 
The   cells   are  numbered   consecutively 
from  the   cell   in  the   northwest   corner 
and   proceeding   from  west   to  east 
southward   (fig.    8). 

2.  Receiving   cell  number:      The   number  of 
the   cell   into  which    the  most   signifi- 
cant  portion  of    the   runoff  drains. 
Drainage   direction  is   determined   by 
cell   topography.      An  example   is   shown 
in  figure  8.      An  arrow  represents   the 
drainage   direction  from  cell  number 
17   as   east,    so   its   receiving   cell 
would  be   number   18.      The   receiving 
cell  number  for   the  watershed   outlet 
must  be   a  number  greater   than  the 
total   number  of   cells   in  the  water- 


6. 


shed.   Thus  the  number  for  the 
receiving  cell  for  cell  number  25  in 
figure  8  would  be  any  number  greater 
than  26. 

SCS  curve  number:   The  runoff  curve 
number  or  hydrologic  soil-cover 
complex  number  used  in  the  SCS  (U.S. 
Department  of  Agriculture,  Soil 
Conservation  Service)  equation  for 
estimating  direct  runoff  from  storm 
rainfall.   Although  the  curve  number 
for  any  antecedent  moisture  condition 
may  be  used,  we  recommend  a  value  for 
moisture  condition  II  to  keep  the 
analyses  consistant.   Table  2  gives 
curve  numbers  for  various  land-use 
conditions  in  rural  areas.   For  other 
conditions  not  listed,  refer  to  the 
Hydrology  Guide  for  Minnesota  (U.S, 
Department  of  Agriculture,  Soil 
Conservation  Service  1976).   For 
moisture  condition  I  or  III, 
equivalent  curve  numbers  are  given  in 
table  3.   If  more  than  one  land-use 
condition  exists  within  a  cell,  use  a 
weighted  average  value  if  possible. 

Land  slope:   The  major  slope,  in 
percent  of  rise,  of  the  cell.   This 
can  be  obtained  from  topographic 
maps,  existing  data,  or  actual  meas- 
urements.  Assume  an  average  slope 
if  the  cell  is  irregular.   If  the 
cell  is  mainly  water  or  marsh,  enter 
a  value  of  0. 

Slope  shape  factor:   An  identifica- 
tion number  used  to  indicate  the 
dominant  slope  shape  of  the  cell: 
1,  a  uniform  slope;  2,  a  convex 
slope;  and  3,  a  concave  slope 
(fig.  10). 

Field  slope  length:   Field  slope 
length  can  be  found  by  using  figure 
11  along  with  table  4.   If  the  cell  is 
mainly  water  or  marsh,  enter  a  value 
of  0. 
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Table   2. 

Runoff    curve   numbers    and    surface-condition  constants   for   various    land-use 

situations 


Surface- 
condition 

Runof  f 

curve 

number 2 

constant-'- 

Soil 

Soil 

Soil 

Soil 

Land-use  condition 

c 

group  A 

group 

B 

group  C 

group  D 

Fallow 

0.22 

77 

86 

91 

94 

Row  crop 

Straight  row 

.05 

67 

78 

85 

89 

Contoured 

.29 

65 

75 

82 

86 

Small  grain 

.29 

63 

74 

82 

85 

Legumes  or  rotation  meadow 

.29 

58 

72 

81 

85 

Pasture^ 

Poor 

.01 

68 

79 

86 

89 

Fair 

.15 

49 

69 

79 

84 

Good 

.22 

39 

61 

74 

80 

Permanent  meadow 

.59 

30 

58 

71 

78 

Woodland 

.29 

36 

60 

73 

79 

Forest  with  heavy  litter 

.59 

25 

55 

70 

77 

Farmsteads 

.01 

59 

74 

82 

86 

Urban  ( 2  1%— 27%  impervious 

.01 

72 

79 

85 

88 

surfaces) 

Grass  waterway 

1.00 

49 

69 

79 

84 

Water 
Marsh 
Animal  lot 

0 
0 

100  ■ 
85  ■ 

Unpaved 

91  ■ 

Paved 
Roof  area 

94  • 
100  ■ 

1 


Source:   Young  et  al .  (1982a). 


-Source:   U.S.  Department  of  Agriculture,  Soil  Conservation  Service  (1976). 
Values  given  are  for  Antecedent  Moisture  Condition  II. 

^Pasture  should  be  considered  "poor"  if  it  is  heavily  grazed  with  no  mulch. 
"Fair"  pasture  has  between  50%  and  75%  plant  cover  and  is  moderately  grazed. 
"Good"  pasture  is  lightly  grazed  and  has  more  than  75%  plant  cover. 
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Table  3. 

Runoff  curve  numbers  for  AMC  I 

and  III  corresponding  to  AMC  II 


Curve 

number  for 
condition 
II 


Curve 

number  for 
conditions — 
I        III 


Table   4. 

Slope    length   in   feet    according    to   major 

field  slope  length  area  and  land  slope 


Major  field 
slope  length 
area 


Slope  length  (feet)  by  slope 
steepness  category  (%) — 
0-2     3-6     7-12     >13 


94 

91 

89 

86 

85 

84 

82 

81 

80 

79 

78 

77 

75 

74 

73 

72 

71 

70 

69 

68 

67 

65 

63 

61 

60 

59 

58 

55 

49 

39 

36 

30 

25 


85 
80 
80 
72 
70 
68 
66 
64 
63 
62 
60 
59 
57 
55 
54 
53 
52 
51 
50 
48 
47 
45 
43 
41 
40 
39 
38 
35 
30 
21 
19 
15 
12 


98 

97 

97 

94 

94 

93 

92 

92 

91 

91 

90 

89 

88 

88 

87 

86 

86 

85 

84 

84 

83 

82 

80 

78 

78 

77 

76 

74 

69 

59 

56 

50 

43 


Source:   U.S.  Department  of 
Agriculture,  Soil  Conservation 
Service  (1972). 


100 

125 

100 

75 

100 

250 

200 

150 

200 

250 

200 

150 

200 

200 

200 

200 

Source:   Otterby  and  Onstad  (1978) 


Uniform 


Convex 


Concave 


Figure  10 

Examples  of  various  slope  shapes  and 
the  identification  numbers  given 
to  each. 


15 


Figure  11 

Division  of  Minnesota  into  major 

field  slope  length  areas. 

Channel  slope:   The  average  slope,  in 
percent,  of  the  channel  or  channels, 
in  the  cell.   This  can  be  obtained 
from  the  topographic  maps,  existing 
data,  or  actual  measurements.   If 
there  is  no  definable  channel  within 
the  cell,  assume  a  series  of  small 
channels,  with  an  average  slope  equal 
to  half  the  land  slope.   If  the  cell 
is  mainly  water  or  marsh,  enter  a 
value  of  0. 

Channel  sideslope:   The  average  side- 
slope,  in  percent,  of  the  channel  or 
channels  in  the  cell.   If  no  value 
can  be  measured  or  estimated,  assume 
a  sideslope  of  10  percent.   If  the 
cell  is  mainly  water  or  marsh,  enter 
a  value  of  0. 

Manning's  roughness  coefficient  for 
the  channel:   Values  of  various 
roughness  coefficients  for  different 
land-use  conditions  at  the  time  of 


the  storm  are  shown  in  table  5.   If 
there  is  no  definable  channel  within 
the  cell,  select  a  roughness  coeffi- 
cient appropriate  for  the  main 
surface  condition  in  the  cell. 

10.  Soil  erodibility  factor:   The  K- 
f actor  used  in  the  USLE,  obtained 
from  SCS  soils  data.   If  the  cell  is 
mainly  water  or  marsh,  enter  a  value 
of  0. 

11.  Cover  and  management:   The  C-f actor 
used  in  the  USLE.   Since  we  are 
working  on  a  storm  basis,  the  value 
to  be  used  is  the  soil  loss  ratio 
corresponding  to  the  appropriate 
period  of  the  growing  season.   In 
order  to  keep  the  analyses  consistent 
among  different  watersheds,  you  will 
normally  examine  the  worst-case 
condition  occurring  during  the  fallow 
or  seedbed  periods.   Appropriate  soil 
loss  ratio  values  can  be  obtained 
from  the  tables  in  appendix  5.   These 
tables  were  copied  in  their  entirety 
from  Wischmeier  and  Smith  (1978).   If 
the  cell  is  mainly  water  or  marsh, 
enter  a  value  of  0.   If  the  cell  is 
mostly  urban  or  residential,  enter  a 
value  of  0.01. 

12.  Support  practice  factor:   The  P- 

f actor  used  in  the  USLE.   To  examine 
a  worst-case  situation,  assume  a 
practice  factor  of  1.0.   Other  values 
for  various  conservation  practices 
can  be  found  in  Wischmeier  and  Smith 
(1978).   P-f actors  for  terraced 
fields  can  be  found  in  table  6.   If 
the  cell  is  mainly  water  or  marsh, 
enter  a  value  of  0.   If  the  cell  is 
mostly  urban  or  residential,  enter  a 
value  of  1.0. 

13.  Surface  condition  constant:   A  value 
based  on  land  use  at  the  time  of  the 
storm  to  make  adjustments  for  the 
time  it  takes  overland  runoff  to 
channelize.   Values  are  shown  in 
table  2. 
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Table  5. 

Manning's  roughness  coefficients  for  channelized  flow 


Natural  channels  *■ 


Description 


Excavated  or  dredged  channels 

Ordinary  concrete  0.013 

Earth,  straight,  uniform,  and  clean  .022 

Same,  but  with  some  short  grass  or  weeds  .027 

Earth,  winding  and  sluggish,  with  no  vegetation  .025 

Same,  but  with  some  grass  or  weeds  .030 

Channels  not  maintained;  weeds  and  some  brush  .080 

Natural  streams 

Clean  and  straight;  no  rifts  or  deep  pools  .030 

Clean  and  winding;  some  pools  and  shoals  .040 

Clean  and  winding;  some  weeds,  stones,  and  pools  .048 

Sluggish  reaches  with  weeds  and  deep  pools  .070 


Cultivated  land  and  waterways 2 


Cover  and  cover  density 


Smooth,  bare  soil 

less  than  1  inch  deep  0.030 

1-2  inches  deep  .033 

2-4  inches  deep  .038 

4-6  inches  deep  .045 

Cornstalks  (assumes  residue  stays  in  place  and  is  not  washed  away) 

1  ton/acre  .050 

2  tons/acre  .075 

3  tons/acre  .100 

4  tons/acre  .130 
Wheat  straw  (assumes  residue  stays  in  place  and  does  not  wash  away) 

1  ton/acre  .060 
1.5  tons/acre  .100 

2  tons/acre  .150 
4  tons/acre  .250 

Grass  (assumes  grass  is  erect  and  as  deep  as  flow) 

Sparse  .040 

Poor  .050 

Fair  .060 

Good  .080 

Excellent  .130 

Dense  .200 

Very  dense  .300 

See  footnotes  at  end  of  table. 
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Table  5. — Continued 

Manning's  roughness  coefficients  for  channelized  flow 


Cultivated  land  and  waterways 2 


Cover  and  cover  density 


Small  grain  (20%  to  full  maturity — rows  with  flow) 

Poor,  7-inch  rows 

Poor,  14-inch  rows 

Good,  7-inch  rows 

Good,  14-inch  rows 
Water  or  marsh3 


0.130 
.130 
.300 
.200 
.990 


1  Source:   Chow  (1959). 

2Source:   Foster  et  al .  (1980). 

^Value  serves  as  a  flag  only  to  tell  the  computer  that  the  surface  is 
water. 


Table  6. 

Sediment  delivery  subf actor,  P,  for 

terraces^ 


Terrace 
grade  (%) 


Closed  outlet^ 
0  (level) 

.1 

,2 

.4 

.6 


>.9 


30.05 
.10 
.13 
.17 
.29 
.49 
.83 
1.00 


Source:   Foster  and  Highfill  (1983). 

Potential  for  net  erosion  in  terrace 

channels  depending  upon  flow  hydraulics 

and  soil  erodibility  in  the  channels.   If 

net  erosion  occurs,  P>1. 

^ Includes  terraces  with  underground 

outlet. 

3Wischmeier  and  Smith  (1978). 


14.  Aspect:   A  single  digit  designating 
the  principal  direction  of  drainage 
from  the  cell.   This  can  be  one  of 
eight  possible  directions,  1  being 
north  and  proceeding  clockwise,  8 
being  northwest  (fig.  9).   If  there 
is  no  drainage  from  the  cell,  input 
a  0. 

15.  Soil  texture:   The  major  soil  texture 
classification  for  the  cell.   The 
major  soil  texture  can  be  determined 
from  the  texture  triangle  shown  in 
figure  12.   The  texture  classes  and 
numbers  to  designate  each  are: 


Texture 

Water 

Sand 

Silt 

Clay 

Peat 


Input 
value 

0 
1 
2 
3 
4 


OaIOO 

0/  \90 

0/  \80 

5>       3  0/  \70 

*?  /  \ 

K        4  0/  >60       ^ 

<?  A  CLAY  -  1.40         /    \  "^ 

^       50/  \  /  \50 

<4  /       \  ' s,|t  cl°y  \         < 

60/    sondy  \- i  \40 

:l°y        ^         clay    loam     \  ="'*    cla* 
\        loam 

'0/     sandy    clay        ^ __"\oL 

loam 
(  loam  /       SILT  -1.30       \20 

SAND   -  1.60 

\  /  ' \ 

,     "^v        sandy    loam   \  J  silt    loam  / 

.  loamy  •^  7  ,      silt 

nd     v  sand    v^^  /  > 

100^ — " =" * *■ " * l ' ^  ° 

0         10       20       30       40       50       60       70       80       90      100 

PERCENT     SILT 

Figure  12 

Soil  bulk  density  (g/cm^)  for 

each  of  the  major  soil  textures. 


16.   Fertilization  level:   A  single  digit 
designation  of  the  level  of  fertili- 
zation on  the  field.   The  number  to 
be  input  and  the  levels  are: 


Assumed 

Level  of 

fertilization 

fertili- 

(lb/acre) 

zation 

N 

P 

Input 

None 

0 

0 

0 

Low 

50 

20 

1 

Medium 

100 

40 

2 

High 

200 

80 

3 

For  a  manure-applied  field,  assume 
low  fertilization  for  an  average 
application  of  manure  and  medium 
fertilization  for  a  heavy  applica- 
tion of  manure.   If  the  cell  is 
mainly  water  or  marsh,  the  fertili- 
zation level  is  0.   If  the  cell  is 
mainly  urban  or  residential,  the 
normal  input  would  be  a  value  of  0. 
But  if  high  levels  of  fertilization 
are  known  to  be  practiced  locally 
(for  example,  lawns),  then  input  an 
appropriate  fertility  level. 


17.  Fertilizer  availability  factor:   The 
percentage  of  fertilizer  left  in  the 
top  half  inch  of  soil  at  the  time  of 
the  storm.   The  worst  case  would  be 
if  none  of  the  fertilizer  had  been 
incorporated  into  the  soil,  an 
availability  factor  of  100  percent. 
Availability  factors  for  various 
tillage  practices  are  shown  in  table 
7.   If  the  cell  is  mainly  water  or 
marsh,  enter  a  value  of  0.   If  it  is 
mostly  urban  or  residential,  enter  a 
value  of  100. 

18.  Point  source  designator:   A  single- 
digit  designator  of  point  sources  in 
the  cell,  such  as  feedlots,  springs, 
and  waste  treatment  plants:   0 
denotes  no  point  sources;  any  other 
number  designates  the  number  of 

Table  7. 

Fertilizer  availability  factors  according 

to  tillage  practice 


Tillage  practices- 


Fertilizer 
availability 
factor  (%) 


Large  offset  disk 

Moldboard  plow 

Lister 

Chisel  plow 

Disk 

Field  cultivator 

Row  cultivator 

Anhydrous  applicator 

Rod  weeder 

Planter 

Smooth 


40 
10 
20 
67 
50 
70 
50 
85 
95 
85 
100 


■'■If  more  than  one  tillage  has  been  made 
since  the  fertilizer  application,  use  the 
product  of  the  two  factors  divided  by 
100. 

Source:   Williams  (1983). 


19 


point  sources  discharging  within  the 
cell  boundaries.   If  one  or  more 
point  sources  are  in  the  cell,  the 
appropriate  point  source  data  will 
need  to  be  input  using  the  format 
specified  in  table  1. 

19.  Gully  source  level:   If  desired,  an 
estimate  can  be  made  for  the  tons  of 
gully  erosion  occurring  in  the  cell. 
This  value  will  then  be  included  in 
the  total  amount  of  sediment  eroded 
in  the  cell. 

20.  Chemical  oyxgen  demand  (COD)  factor: 
A  value  for  the  COD  concentration 
from  the  cell,  based  on  the  land  use 
in  the  cell.   The  COD  factors  for 
various  land-use  situations  are 
shown  in  table  8. 

21.  Impoundment  factor:   A  factor 
indicating  the  presence  of  an 
impoundment  terrace  system  within 
the  cell.   Zero  would  indicate  no 
terrace  in  the  cell;  any  other 
number  would  be  the  number  of 
impoundments  in  the  terrace  system, 
with  a  maximum  of  13.   The  area  in 
acres  draining  into  each  impoundment 
and  the  diameter  in  inches  of  the 
outlet  pipe  of  each  impoundment  will 
need  to  be  input  using  the  format 
specified  in  table  1. 

22.  Channel  indicator:   A  single  digit 
indicating  the  presence  of  a  defined 
channel  within  the  cell:   0  denotes 
no  defined  channel;  any  other  number 
indicates  the  number  of  channels  in 
the  cell. 


Point  Source  Inputs 

Nonfeedlot:   If  there  is  a  point  source 
in  a  cell  that  is  not  a  feedlot,  such  as  a 
spring  or  a  waste  treatment  plant  dis- 
charge, the  information  describing  the 
source  must  come  immediately  after  the 
basic  cell  data.   As  described  in  table  1, 


the  input  information  includes  a  1  as  the 
first  digit,  to  indicate  that  the  source 
is  not  a  feedlot.   The  remaining  informa- 
tion includes  the  incoming  flow  rate  in 
cubic  feet  per  second  and  the  inflow 
concentration  of  N,  P,  and  COD  in  parts 
per  million.   These  data  are  entered  as  a 
separate  line  for  each  nonfeedlot  point 
source  in  the  cell. 

Feedlots:   In  the  event  that  the  point 
source  input  in  any  of  the  cells  in  the 
watershed  is  from  an  animal  feedlot,  the 
feedlot  variables  mentioned  in  the  data 
file  format  (table  1)  will  have  to  be 
input.   The  first  digit  will  be  a  2  to 
indicate  that  the  source  is  a  feedlot. 
For  an  explanation  of  the  remaining 
variables,  refer  to  Young  et  al .  (1982a). 
All  variables  required  for  input,  except 
the  nitrogen  factor,  can  be  found  in  Young 

Table  8. 

Chemical  oxygen  demand  (COD)  factors  for 

various  land-use  situations 


Land  use 


1 


COD  factor  (mg/L) 


Row  crops 

Small  grain 

Pasture  and  open 

Alfalfa 

Forested 

Fallow 

Farmsteads  and  urban 

nonresidential 
Water 
Marsh 


170 

80 
60 
20 
65 
115 
80 

0 
25 


1  Sources  of  data  are  as  follows:   Row 
crops  and  fallow,  Thompson  et  al .  (1978), 
Harms  et  al .  (1974);  small  grain  and 
alfalfa,  Harms  et  al.  (1974);  pasture  and 
open  land,  Crow  et  al .  (1979),  Thompson 
et  al.  (1978),  Harms  et  al.  (1974); 
forested  land,  Timmons  et  al .  (1977),  R. 
A.  Young,  unpublished  data;  and  farmstead 
and  urban  nonresidential,  Weibel  (1969). 
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et  al.    (1982a).      The   nitrogen   factors   can 
be   obtained   from  table   9.      The   feedlot 
data  can  be   recorded   as   shown   in  figure  3. 


Impoundment  Terraces 

Item  number  21  indicates  the  number  of 
impoundments  in  the  cell.   The  data 
describing  the  impoundments  will  come 
immediately  after  the  basic  cell  data, 
described  in  table  1,  the  impoundment 
variables  include  the  area  draining  into 
the  impoundment  in  acres  and  the  diameter 
of  the  impoundment  outlet  pipe  in  inches. 
These  variables  are  entered  in  pairs, 
describing  each  impoundment  in  the  terrace 
system.   The  impoundment  data  can  be 
collected  as  shown  in  figure  4. 


As 


Table  9. 

Ratio  of  total  nitrogen  produced  by 
various  animals  to  that  produced  by  a 
1,000-pound  slaughter  steer 


Design 

weight^ 

Ratio 

Animal  type^ 

(pounds) 

N 

Slaughter  steer 

1,000 

1.00 

Young  beef 

500 

.60 

Dairy  cow 

1,400 

1.68 

Young  dairy  stock 

500 

.46 

Swine 

200 

.26 

Feeder  pig 

50 

.07 

Sheep 

100 

.13 

Turkey 

10 

.02 

Chicken 

4 

.01 

Duck 

4 

.01 

Horse 

1,000 

.81 

iData  from  Midwest  Plan  Service  (1975) 
except  for  swine,  which  is  from  American 
Society  of  Agricultural  Engineers  (1982). 

^Interpolation  of  values  should  be  based 
on  the  maximum  weight  animals  would  be 
expected  to  reach. 


Recording  Data 

The  data  file  can  be  collected  and 
organized  on  data  sheets.   Separate 
sheets  are  provided  for  collection  of  the 
watershed  data  (fig.  2),  the  cell  data 
(fig.  1),  the  impoundment  data  (fig.  4), 
and  the  feedlot  data  (fig.  3).   In  most 
cases,  general  information  about  the  soil 
characteristics  and  the  land  use  in  each 
cell  must  be  collected  before  many  of  the 
cell  parameters  can  be  established. 
Supplementary  data  collection  sheets  are 
provided  for  this  purpose  (fig.  5).   It 
should  be  noted  that  values  given  in  the 
tables  and  figures  are  based  mainly  on 
average  values  found  in  the  literature. 
If,  for  any  variable,  you  recognize  that 
the  value  given  in  the  table  or  figure  is 
not  correct  for  that  particular  cell, 
input  a  value  you  know  to  be  more 
accurate. 

For  those  cells  where  you  cannot  decide 
what  value  to  input  for  a  particular 
variable  because  of  nonuniform  conditions 
in  a  cell,  use  a  weighted  average 
whenever  possible.   If  this  is  not 
practical  because  of  time  or  labor 
limitations,  a  value  should  be  selected 
corresponding  to  the  predominant 
condition  in  the  cell.   This  applies  to 
parameters  3,  4,  7,  8,  9,  10,  11,  12,  13, 
17,  and  20. 

Appendix  2  contains  a  demonstration  of 
how  to  set  up  a  data  file  for  input  into 
the  model. 


Using  MLMIS  Data 

A  data  file  containing  many  of  the  cell 
parameters  required  as  input  can  be 
obtained  from  MLMIS  (see  "Introduction  to 
User's  Guide").   Using  the  MLMIS  data 
file^  can  minimize  the  actual  number  of 


^An  explanation  of  the  MLMIS  data  file 
can  be  found  in  appendix  3. 
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parameters  you  must  search  for.   Collect 
the  information  not  provided  by  MLMIS  and 
record  it  in  a  computer  file  called  an 
edit  file.   Along  with  these  data,  you  may 
also  enter  values  into  the  edit  file  to 
override  the  MLMIS  data.   A  separate 
computer  program  (DBDFL)  has  been  written 
to  combine  the  two  files  and  form  the 
required  model  input  data  file.   The 
requirements  for  establishing  a  data  file 
through  this  method  are: 

1.   A  data  file  on  the  watershed  obtained 
from  MLMIS,  which  contains  the 
following  information: 

a.  row  designation 

b.  column  designation 

c.  township  number 

d.  range  number 

e.  section  number 

f.  40-acre  parcel/or  10-acre  parcel 

g.  minor  watershed 
h.  primary  land  use 
i.  forest  cover 

j.  soil  erodibility  code 

k.  land  slope  (percent) 

1.  soil  texture 

m.  river  slope  (percent) 

n.  30°  aspect 

o.  river-channel  designation 

p.  eight-way  aspect 

The  base  level  of  information  used  in 
this  data  file  must  be  compatible 
with  the  one  you  set  up.   That  is, 
the  cell  division  established  must  be 
identical  for  both  data  sources.   The 
MLMIS  data  file  begins  in  the  north- 
west corner  of  the  watershed  and 
proceeds  in  the  same  fashion  as  the 
numbering  system  described  in  step  4 
of  the  preliminary  examination. 


4.  A  desired  area  (acres)  for  each  cell. 

5.  Storm  characteristics.  The  user  has 
the  option  of  using  a  default  storm 
such  as  one  of  25-year,  24-hour 
frequency,  or  of  entering  different 
storm  characteristics.  If  the  user 
chooses  to  enter  a  different  storm, 
the  following  data  will  be  required: 

a.  storm  precipitation  (inches)^ 

b.  storm  energy-intensity  valued 

6.  Field  slope  length  area  number.   This 
number  can  be  obtained  from  figure  11. 

7.  If  no  value  has  been  recorded  in  the 
MLMIS  data  file  for  the  soil  erodi- 
bility of  the  cell,  one  will  have  to 
be  input. 

In  using  DBDFL  to  generate  a  parameter 
file,  you  can  input  as  many  known  values 
as  desired.   DBDFL  takes  the  information 
previously  mentioned  and  generates  a 
preliminary  data  file.   This  contains 
variable  estimates,  or  default  values, 
made  by  the  computer  using  the  available 
information.   Since  the  computer  can  only 
make  rough  estimates  of  the  variables,  we 
do  not  recommend  watershed  analysis  using 
the  default  file.   Through  use  of  a 
separate  file  containing  parameter  values 
you  have  chosen,  you  can  edit  the 
preliminary  data  file.   Through  the 
editing  process,  you  can  input  as  many 
known  values  into  the  final  data  file  as 
desired. 

A  listing  of  the  parameters  in  the  edit 
file  is  shown  in  table  10.   Item  3, 
land-use  code,  can  be  determined  using 
table  11.   All  other  items  are  as 


2.  A  name  to  be  applied  to  the  output 
data  file  (up  to  six  characters). 

3.  A  30-character  name  to  be  applied  to 
the  watershed  and  the  data  file. 


^Refer   to   appendix    1    for   rainfall 
frequency   and   energy-intensity   curves, 


2  2 


Table  10. 

Parameter  explanation  for  the  edit  file 


Edit  file 


Line  1  to  end  of  file — Cell  parameters 

1.  Cell  number 

2.  Aspect^ — number   indicating   direction  of   drainage   (fig.    9) 

3.  Land-use   code    (table    11) 

4.  SCS   curve   number   (table    2) 

5.  Land   slope    (%) 

6.  Slope  shape  factor  (fig.  10) 

1 — uniform  slope 
2--convex  slope 
3 — concave  slope 

7.  Channel   slope    (%) 

8.  Channel    sideslope    (%) 

9.  Manning's   roughness  coefficient    for  the   channel    (table   5) 

10.  Soil   erodibility   factor 

11.  Cropping   factor 

12.  Practice   factor 

13.  Fertilization  level 

1 — low  fertilization 
2 — average  fertilization 
3 — high  fertilization 
4 — no  fertilization 

14.  Fertilizer  availability  factor  (%)  (table  7) 

15.  Point  source  designator  (a  single  digit  number  other  than  0 
designates  the  number  of  sources  in  the  cell) 

16.  Gully  source  level  (tons) 

17.  Impoundment  factor  (indicates  the  number  of  impoundments) 

Format  (14,12, 13, 14, F5. 1 , 12 ,2F5. 1 ,F5.3 ,3F4.2 , 12 , 14 ,12, 14, 13) 

*At  the  present  time,  the  aspects  given  by  the  MLMIS  data  file  have  been 
found  to  be  inaccurate.   For  this  reason,  input  of  the  aspects  of  all 
of  the  cells  into  the  edit  file  is  required. 

described  in  "Establishing  a  Watershed  DBDFL  and  the  various  options  available, 

Data  File."   In  the  edit  file,  you  may  including  the  edit  file,  are  shown  in 

input  the  values  for  whichever  items  are  appendix  4. 

known  and  leave  the  rest  blank.   When  a  0 

value  or  a  blank  is  found  by  the  computer,  A  collection  form  for  the  edit  file  can 

no  change  will  be  made  for  that  item  in  be  found  in  figure  6. 

the  data  file.   An  example  of  the  use  of 
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Table  11 

Land-use  codes  for  the  edit  file, 


Land-use  category 


Fallow 
Row  crop 

Straight  row 

Contoured 
Small  grain 

Legumes  or  rotation  meadow 
Pasture 

Poor 

Fair 

Good 
Permanent  meadow 
Woodland 

Forest  with  heavy  litter 
Farmsteads 


Code 


2 
3 
4 
5 

6 

7 

8 

9 

10 

11 

12 


Model  Output 

The  basic  output  from  AGNPS  includes 
hydrology,  with  estimates  of  both  volume 
and  peak  runoff;  sediment,  with  estimates 
of  upland  erosion,  channel  erosion,  and 
sediment  yield;  and  nutrients,  with 
estimates  for  nitrogen  (N),  phosphorus 
(P),  and  chemical  oxygen  demand  (COD)  in 
concentration  and  mass  units.   Sediment 
results  are  available  for  five  particle- 
size  classes:   clay,  silt,  small  aggre- 
gates, large  aggregates,  and  sand. 

AGNPS  allows  various  output  options.   The 
outputs  can  be  examined  for  a  single  cell 
or  for  the  entire  watershed.   A  prelimi- 
nary output,  given  for  all  watersheds 
being  examined,  includes  watershed  and 
cell  areas,  storm  precipitation  and  ero- 
sivity,  estimated  values  at  the  watershed 
outlet  of  runoff  volume  and  peak  flow 
rate,  and  a  detailed  analysis  of  the 
sediment  and  nutrient  yields  (fig.  13). 
The  detailed  sediment  analysis  includes 
area-weighted  erosion  rates  for  both 
upland  and  channel,  sediment  delivery 


ratios,  sediment  enrichment  ratios,  mean 
sediment  concentrations,  area-weighted 
yields,  and  net  sediment  yields.   These 
values  are  given  for  each  of  the  five 
particle-size  classes  of  sediment,  as  well 
as  for  the  total  sediment.   The  detailed 
nutrient  analysis  includes  the  N  and  P 
mass  per  unit  area  for  sediment-adsorbed 
nutrients;  the  soluble  N,  P,  and  COD  mass 
per  unit  area  in  runoff;  and  the  N,  P,  and 
COD  concentration  in  the  runoff. 

Additional  information  is  also  available 
for  selected  individual  cells  in  the 
watershed  when  requested  (fig.  14). 
Information  given  when  examining  individ- 
ual cells  includes  runoff  and  sediment 
analyses  for  each  cell.   Runoff  analysis 
provides  estimates  for  each  cell  of  drain- 
age area,  runoff  volume,  percent  of  runoff 
volume  entering  the  cell  from  above,  and 
peak  runoff  rate.   Sediment  analysis  pro- 
vides estimates  for  each  cell  of  upland 
erosion  rate,  amount  of  sediment  generated 
within  the  cell,  amount  entering  the  cell 
from  above,  the  sediment  yield  leaving  the 
cell,  and  percent  deposition  in  the  cell. 
A  negative  sign  preceding  percent  of 
deposition  indicates  channel  scouring  in 
the  cell.   An  option  also  exists  to 
provide  the  cell  sediment  analysis  broken 
down  into  the  five  particle-size  classes: 
clay,  silt,  small  aggregates,  large 
aggregates,  and  sand. 

A  detailed  nutrient  analysis  for  each 
cell  is  also  available  that  provides 
estimates  of  adsorbed  and  soluble 
nutrients  in  mass  per  unit  area  and  the 
concentration  of  these  nutrients  in  the 
runoff. 

Execution  of  AGNPS  simply  requires  the 
input  of  a  data  file  name  and  choosing 
the  output  options  desired.   Detailed 
sediment  and  nutrient  options  provide  a 
great  deal  of  information  and  should  only 
be  requested  if  a  very  extensive  investi- 
gation of  the  watershed  is  taking  place. 
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*  AGNPS  Version  2.2       Input  file:  MB045 

*  AGRICULTURAL  NON- POINT  SOURCE  POLLUTION  MODEL 

*  Watershed  studied:  MEADOW  BROOK  WTSHD#045  1983 

x  x  x  x  x  it  x  x  x  ^  x  x  x  x  x  x  x  x  x  x  x  x  -A'  ■£■  x  x  x  x  x  x  x  x  x  x  x  x  x  x  x  x  x  x  x  x  vV  x  x  x  x  x  x  x  x  >>  "A"  'A-  vV  x  x  x  x  'A-  x  x  x  x  'A  -A-  x  x  x  "A-  -A-  -A-  ■ 

The  area  of  the  watershed  is:  8680  acres 

The  area  of  each  cell  is:  40.0  acres 

The  characteristic  storm  precipitation  is:  4.4  inches 

The  storm  energy- intensity  value  is:  56 
The  cell  parameters  are  derived  for: 

OUTFLOW  FROM  FEEDLOT    1 
Nitrogen  concentration  at  the  discharge  point  (ppm)    112.2 
Phosphorous  concentration  at  discharge  point  (ppm)     55.12 

COD  concentration  at  discharge  point  (ppm)  1622. 

Nitrogen  mass  at  the  discharge  point  (lbs)  298.7 

Phosphorous  mass  at  discharge  point  (lbs)  146.7 

COD  mass  at  discharge  point  (lbs)  4318. 


==  VALUES  AT  THE  WATERSHED  OUTLET,  CELL  NUMBER   213 


RUNOFF: 

Runoff  volume  (in.): 
Peak  runoff  rate  (cfs): 

NUTRIENTS: 


2.37 
3916 


Total  nitrogen  in  sediment  (lbs/acre) : 

Total  soluble  nitrogen  in  runoff  (lbs/acre): 

Soluble  nitrogen  concentration  in  runoff  (ppm) : 

Total  phosphorus  in  sediment  (lbs/acre): 

Total  soluble  phosphorus  in  runoff  (lbs/acre): 

Soluble  phosphorus  concentration  in  runoff  (ppm) : 

Total  soluble  chemical  oxygen  demand  (lbs/acre): 

Soluble  chemical  oxygen  demand  concentration  in  runoff  (ppm) 

SEDIMENT  ANALYSIS: 


3.16 

1.92 

3.6 

1.58 

.37 

.7 

66.24 

123 


Area  Weighted  Area 

Erosion     Delivery  Mean  Weighted 

Particle  Upland  Channel  Ratio  Enrichment  Concentration  Yield    Yield 

Type     (t/a)  (t/a)     (%)  Ratio         (ppm)  (t/a)    (tons) 


CLAY 

.17 

.00 

99 

3.0 

SILT 

.26 

.00 

80 

3.0 

SAGG 

1.62 

0.00 

38 

1.0 

LAGG 

1.00 

.00 

0 

0.0 

SAND 

.19 

.00 

1 

0.0 

TOTL 

3.23 

0.00 

31 

1.0 

615 

.2 

1433.8 

769 

.2 

1793.5 

2319 

.6 

5406.6 

12 

.0 

27.5 

4 

.0 

8.6 

3719 

1.0 

8670.0 

Figure  13 

AGNPS  output  for  a  8,680-acre  watershed. 
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MODEL  DOCUMENTATION 


Introduction  to  Model  Documentation 

This  section  provides  the  technical  back- 
ground for  the  water  quality  model  AGNPS. 
AGNPS  was  written  in  Fortran  IV  and  devel- 
oped on  a  Hewlett-Packard  1000  computer 
system.  The  program  developed  on  that 
system  was  then  modified  to  run  on  a 
Hewlett-Packard  41  CV  hand-held  calculator 
system  for  use  with  smaller  watersheds. 

This  section  also  discusses  the  computer 
program  DBDFL,  which  is  operated  in  con- 
junction with  AGNPS.   DBDFL  was  devel- 
oped to  assist  the  user  in  creating  model 
input  data  files  for  the  computer  system. 

The    calculations  made   by  AGNPS   occur    in 
three  stages,    or  loops.      Initial   calcula- 
tions  for  all   cells    in  the  watershed   are 
made    in   the   first    loop.      These    calcula- 
tions   include   estimates    for  upland 
erosion,   overland   runoff  volume,    time 
until    overland   flow  becomes    concentrated, 
level   of    soluble   pollutants   leaving   the 
watershed   via  overland    runoff,    sediment 
and    runoff   leaving   impoundment-terrace 
systems,    and   pollutants   coming   from  point- 
source   inputs   such   as    tile    lines    or 
f eedlots. 

The   second   loop   calculates    the   runoff 
volume   leaving   the   cells   containing 
impoundments   and    the   sediment   yields    for 
primary   cells.      A  primary   cell   is   one   that 
no   other   cell   drains    into.      The   sediment 
from  these   and   other  cells   is   broken  down 
into   five   particle-size   classes:      clay, 
silt,    small   aggregates,    large    aggregates, 
and   sand.      Table    12   shows    the   diameter 
range   for  each  particle   size,    the  densi- 
ties,   the  particle   fall  velocities,    and 
the   equivalent    sand   diameters. 

The  sediment  and  nutrients  are  routed 
through  the  rest  of  the  watershed  in  loop 
3.   Calculations  are  made  to  establish  the 
concentrated  flow  rates,  to  derive  the 
channel  transport  capacity,  and  to  calcu- 
late the  actual  sediment  and  nutrient  flow 
rates . 


Preliminary  Calculations 

Base  values  for  each  of  the  cells  in  the 
watershed  are  established  in  loop  1  of 
the  programs.   Estimates  are  made  to 
determine  overland  runoff  volume  and 
duration,  erosion,  and  soluble  pollutant 
flow  occurring  in  each  of  the  cells.   The 
first  execution  made  in  this  loop  is  to 
read  in  the  data  file  for  the  watershed. 
The  parameters  of  the  data  files  are 
shown  in  table  1. 

Total  detached  sediment  and  the  breakdown 
into  the  particle-size  classes  are  made 
on  a  per  cell  basis.   A  modified  USLE  is 
used  to  predict  upland  erosion  for  single 
storm  events  (Wischmeier  and  Smith 
1978). 


E=EI*Ks*Lf*Sf*Cf*Pf*SSF, 
where 


[1 


E=soil    loss   in  tons/acre, 
EI=rainfall  energy-intensity   in  hundred 

foot-ton  inch/acre  hour, 
Ks=soil   erodibility   factor   in  ton-acre 

hour/hundred-acre   foot-ton  inches, 
Lf=slope-length  factor, 
Sf =slope-steepness   factor, 
Cf=cover  and  management   factor, 
Pf=support   practice   factor,   and 
SSF=a   calculated   factor  to   adjust   for 

slope   shape. 

The  variables  EI,  Ks ,  Lf ,  Sf ,  Cf ,  and  Pf 
are  as  defined  in  Wischmeier  and  Smith 
(1978).   The  slope  shape  factors  were 
calculated  using  complex  slope  factors 
found  in  Wischmeier  and  Smith  (1978).   The 
factors  were  developed  to  adjust  the 
erosion  estimates  for  convex  and  concave 
slopes.   Calculations  for  a  convex  slope 
were  based  on  a  75-foot  slope  the  upper 
third  of  which  had  a  gradient  of  2  per- 
cent, the  middle  third  7  percent,  and  the 
lower  third  12  percent.   For  the  concave 
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slope,  a  75-foot  slope  was  again  used, 
with  the  upper  third  having  a  gradient  of 
12  percent,  the  middle  third  7  percent, 
and  the  lower  third  2  percent.   Table  13 
shows  the  slope  shapes  and  the  factors 
used  for  each. 

Soil  texture  is  used  to  determine  the 
fractional  distribution  of  the  eroded 
sediment  into  particle  classes  (table  14). 
The  runoff  volume  from  each  cell  is 
determined  using  the  U.S.  Department  of 
Agriculture,  Soil  Conservation  Service 
(1972)  curve  number  method.   The  equation 
is 


RF= 


.(RL-0.2*S)2 
RL+0.8*S 


[2 


where 

RF=runoff   in  inches, 

RL=storm  precipitation  in  inches,  and 
S=a  retention  factor  in  inches. 

The  retention  factor  is  defined  as 

S=U,000/CN-10,  [3; 

where 

CN  is  the  curve  number. 

The  curve  number  is  a  cell  characteristic 
dependent  on  land  use,  soil  type,  and 
hydrologic  soil  conditions. 

The  overland  runoff  duration,  or  the  time 
needed  for  concentrated  flow  to  occur,  is 
calculated  using  the  runoff  velocity  as 
determined  in  U.S.  Department  of 
Agriculture,  Soil  Conservation  Service 
(1972). 

The  equation  used  for  computation  is 

0FT=Ls/Vo,  [4] 


where 

OFT=overland  flow  time  in  seconds, 
Ls=field  slope  length  in  feet,  and 
V0=overland  flow  velocity  in  feet  per 
second . 

The  velocity  is  calculated  as 
0.5*LOG   (S  *100)-SCC 

v0=io     io  i 

where 


[5: 


S;L=land    slope    in  feet/foot,    and 
SCC=overland   surface    condition  constant. 

The  surface  condition  constant  is  a  cell 
characteristic  that  accounts  for  the 
effects  of  land  use  and  vegetation. 

The  methods  used  to  predict  N  and  P  yields 
from  the  cells  and  watershed  were 
developed  by  Frere  et  al .  (1980)  and 
appear  in  CREAMS  (Chemicals,  Runoff,  and 
Erosion  from  Agricultural  Management 
Systems).   The  contributions  of  soluble  N 
and  P  from  each  of  the  cells  are 
calculated  in  loop  1  and  routed  into  the 
channel.   Once  the  soluble  nutrients  reach 
concentrated  flow,  they  are  assumed  to 
remain.   That  is,  the  amount  arriving  in 
the  overland  flow  from  any  particular  cell 
is  simply  added  to  what  is  already  present 
in  the  channel,  with  no  losses  of  soluble 
nutrients  in  the  channel  allowed. 

Figure  15  illustrates  the  concepts  used  in 
the  nutrient  portion  of  the  model  with 
regard  to  soluble  forms  of  N  and  P  in 
runoff  waters.   As  the  figure  shows,  the 
soluble  N  and  P  calculations  account  for 
the  effects  of  rainfall,  fertilization, 
and  leaching.   The  basic  equation  to 
predict  the  soluble  N  and  P  concentration 
in  the  runoff  is 


R0-=C*EXK"2*RO*0.01, 


[6] 
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Table    12. 

Particle-size   breakdown   and   physical   description 


Equivalent 
Particle  Fall  sand 

range  Density  velocity  diameters 

Particle  (mm)  (g/cm3)  (ft/s)  (micrometers) 


Clay  <0.002  2.60  0.102  *   10-4  2 

Silt  0.002-.050  2.65  .263  *    10"3  10 

Small  aggregates  .020-. 075  1.80  .125  *    10~2  20 

Large   aggregates  .200-1.000  1.60  .542  *   10"1  158 

Sand  .050-2.000  2.65  .759   *    10"1  201 


Source:      Foster  et   al.    (1980). 

Table  13. 

Slope  shapes  and  calculated  universal  soil 

loss  equation  multiplication  factors 


Slope  type  SSF 


Uniform  1.00 

Convex  1.30 

Concave  .88 


•^Slope  shape  multiplication  factor. 

Table  14. 

Sediment  particle-size  distribution  (%)  of  the  major  soil  texture  groups 


Small         Large 
Soil  texture       Clay      Silt    aggregates    aggregates       Sand 


High  clay  soils  0.10 

Silt  soils  .05 

High  sand  soils  .02 

Peat  soils  1.00 


Source:   Foster  et  al.  (1980). 
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mg    nut 

ri- 

where 

RO~=N  or  P  concentration  in  the 
runoff , 
C=mean  concentration  of  the  soluble 
portion  of  the  nutrient  in  the  soil 
surface  during  runoff, 
EXK~2=an  extraction  coefficient  for 
movement  into  the  runoff;  and 
RO=total  runoff  in  millimeters. 


The  equation  used  to  predict  soluble  N 
concentration  in  the  model  is 

RON=0.892*((CZERON-CHECKN)*EXP(-XKFN1   [7; 
*EFI)-(CZER0N-CHECKN)*EXP(-XKFN1 
*EFI-XKFN2*R0))/C0EFF+RN*R0/EFRAIN, 

where 

RON=soluble   N   in   the    runoff    in 
pounds   per  acre, 
CZERON=available    soluble   N   content    in 

the   soil   in  kilograms   per  hectare, 
CHECKN=available   N  due    to    the    rainfall   in 
kilograms    per  hectare, 
XKFNl=a   rate   constant    for  downward   move- 
ment  of   N  into   the   soil, 


EFI=total    infiltration   for   the   storm 
in  millimeters, 
XKFN2=a   rate   constant    for  N  movement    into 
the   runoff 
RO=total    storm   runoff   in  millimeters, 
COEFF=a   porosity   factor, 

RN=N   contribution  due    to    tVie    rain   in 
kilograms    per  hectare,    and 
EFRAIN=ef fective   rainfall   in  millimeters. 

The   available   N  in  the   soil   is   calculated 
by   accounting   for   organic   matter  N, 
fertilizer  N,   and   soil   porosity: 


CZERON=(SOLN+FN(X)*FA(X))*COEFF, 
where 


[8] 


SOLN=soluble  N  in  the  surface  centimeter 

of  the  original  soil  in  kilograms 

per  hectare, 
FN(X)=N  fertilizer  application  in  cell  X 

in  kilograms  per  hectare,  and 
FA(X)=fraction  of  this  application 

remaining  in  the  top  centimeter  of 

the  soil. 

The  initial  soluble  N  in  the  top  centi- 
meter of  the  original  soil  is  estimated  by 


SOLN=0.10*CSN*POR, 


[9] 


where 


CSN=concentration  of  N  in  the  pore  water 
of  the  surface  centimeter  of  soil  in 
parts  per  million,  and 

POR=soil  porosity. 

Because  very  few  areas  have  the  soluble 
nutrient  values  for  the  original  soils,  a 
CSN  value  of  5  p/m  from  an  observed  range 
of  2  p/m  to  5  p/m  (Frere  et  al .  1980)  is 
assumed.   The  soil  porosity  is  determined 
by 


P0R=l-(bulk  density/2.65). 


10 
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The   soil    bulk  density,    in  grams    per   cubic 
centimeter,    is   determined   using    the   major 
soil    texture   of    the   cell   as    shown   in 
figure    12.      The    fertilizer   application 
and    the   fraction   remaining   are   both 
variable   inputs.      The   porosity  factor, 
COEFF,    is   solved   by  using   the  porosity  as 

C0EFFO.00001/P0R.  [IT 

The   available   N  due   to   the   rainfall   is 
solved   by   using    the    relation 

CHECKN=RCN*1.0E-06,  [12; 

where 

RCN=N   concentration  in   the   rainfall. 

A  value  of  0.8  mg/L  was  chosen  for  use  in 
the  model. 

The  equation  for  the  rate  constant  for 
downward  movement  of  soluble  N  is 


XKFN1=EXKN1/(10*POR), 
where 


[13] 


EXKNl=extraction  coefficient  for  downward 
movement . 

EXKN1  is  assumed  constant  and  equal  to 
0.25. 

The  effective  infiltration  is  defined  as 

EFI=EFRAIN-RO,  [14] 

where 

EFRAIN=ef fective   rainfall   in  millimeters, 
and 
R0=total   storm  runoff   in  millimeters. 

The   EFRAIN  term  is    calculated   using    the 
equation 


where 

R=storm  rainfall  in  millimeters. 

The  rate  constant  for  movement  of  soluble 
N  into  the  runoff  is 

XKFN2=EXKN2/(10*P0R),  [16] 

where 

EXKN2=extraction  coefficient  for 
movement  into  runoff. 

EXKN2  is  assumed  constant  and  equal  to 
0.050. 

The  N  contribution  due  to  the  rain,  RN, 
is  calculated  as 


RN=RCN*R*0.01. 


[17] 


The   equation  used    to   predict    soluble    P   in 
the   runoff   is   similar  to   the   equation  for 
N   except    that   the   effects   of    rainfall   are 
omitted.      Since   very  little   soluble   P   is 
found    in   rainfall,    its   net    contribution   to 
the    runoff    is   negligible.      The   equation   is 

ROP=0.892*((CZEROP-CHECKP)*EXP(-XKFP1 
*EFI)-(CZER0P-CHECKP)*EXP(-XKFP1 
*EFI-XKFP2*RO))/COEFF+CHECKP*XKFP2 
*R0/C0EFF),  [18] 


where 

ROP: 
CZEROP; 

CHECKS 
XKFPL 
XKFP2: 


=soluble  P  in  the  runoff  in  pounds 
per  acre, 

=available  P  due  to  natural  and 
fertilizer  nutrient  level  in 
kilograms  per  hectare, 

=available  P  due  to  the  initial 
soil  in  kilograms  per  hectare, 

-a   rate  constant  for  downward  move- 
ment of  P  into  the  soil,  and 

=the  rate  constant  for  movement 
into  runoff. 


EFRAIN=R-(10*P0R), 


[15] 
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The  available  soluble  P  in  the  soil  is 
calculated  as 

CZEROP=(SOLP+FP(X)*FA(X) )*COEFF,     [ 19] 

where 

SOLP=soluble   P  in  the   top   centimeter  of 
the  original   soil   in  kilograms   per 
hectare,   and 
FP(X)=P  fertilizer  application  in  cell  X 
in  kilograms   per  hectare. 

The  initial  soluble  P  in  the  top  centi- 
meter of  the  original  soil  is  solved  by 
using   the   relation 


SOLP=0.10*CSP*POR, 
where 


20] 


CSP=concentration  of  P  in  the  pore  water 
of  the  surface  centimeter  of  soil  in 
parts  per  million. 

A  CSP  value  of  2  p/m  is  assumed  in  the 
model.   The  COEFF  term  is  the  same  as 
described  for  the  soluble  N  equation. 

The  available  P  due  to  the  initial  soil 
is  solved  as 

CHECKP=SOLP*COEFF.  [21] 

The  equation  for  the  rate  constant  for 
downward  movement  of  P  is 

XKFP1=EXKP1/(10*POR),  [22; 

where 

EXKPl=extraction  coefficient  for  downward 
movement,  assumed  equal  to  0.25. 

The  rate  constant  for  movement  into  the 
runoff  is 


where 

EXKP2=extraction  coefficient  for  movement 
into  runoff. 

EXKP2  is  also  assumed  constant  and  equal 
to  0.025. 

Chemical  oxygen  demand  (COD)  in  the 
model  is  assumed  soluble.   Calculations 
of  the  amount  of  soluble  COD  in  the 
runoff  are  based  on  the  runoff  volume  and 
the  average  concentration  of  COD  in  that 
volume.   Various  background  concentra- 
tions of  COD  of  runoff  waters  are  shown 
in  table  15.   Estimates  made  from  these 
values  are  used  as  a  basis  for  predicting 
COD  concentration  in  runoff  from  each 
cell.   Soluble  COD  is  assumed  to 
accumulate  only  after  flow  has  become 
channelized,  without  any  allowable 
losses. 

Sediment  and  runoff  routing  through 
impoundment  terrace  systems  is  simulated 
using  relationships  developed  by  Laflen 
et  al.  (1978)  and  Foster  et  al.  (1980). 
These  relationships  were  developed  for 
terrace  systems  that  concentrate  the 
overland  flow  into  pipe-outlet  ponds. 
The  fraction  of  each  specific  particle 
class  that  passes  through  the  impoundment 
is  expressed  as 


Fpi=A1[EXP(B1*Du)-EXP(B1*Ds)]/ 

(B^AD), 

where 

Fpi=fraction  of    that   particle   class 
passing   through    the   impoundment, 
Aj=a   coefficient, 
B^=a   coefficient, 

Du=equivalent  sand  diameter  for  that 
particle  class  in  micrometers, 


[24; 


XKFP2=EXKP2( 10*POR) , 


[23] 
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Table  15. 

Background  concentrations  of  chemical  oxygen  demand  of  runoff  waters 

from  different  land  sources 


Source 


COD 
(mg/L) 


Native  prairie — Minnesota 
Grassland — Michigan 
Grazed  prairie — Oklahoma 
Aspen-birch  forest — Minnesota 
Aspen-birch  forest — Minnesota 
Corn  stubble — Michigan 
Corn — South  Dakota 
Fallow — Michigan 
Fallow — South  Dakota 
Alfalfa — South  Dakota 
Pasture — South  Dakota 
Oats — South  Dakota 
Urban — Ohio 


48 

98 

36 

54 

78 

121 

225 

118 

109 

22 

49 

82 

111 


Reference 


Timmons  and  Holt  (1977) 
Thompson  et  al.  (1978) 
Crow  et  al.  (1979) 
Timmons  et  al.  (1977) 
Young,  unpublished  data 
Thompson  et  al.  (1978) 
Harms  et  al.  (1974) 
Thompson  et  al.  (1978) 
Harms  et  al.  (1974) 
Harms  et  al .  (1974) 
Harms  et  al.  (1974) 
Harms  et  al.  (1974) 
Weibel  (1969) 


Ds=equivalent    sand   diameter   for  the   next 
smallest   particle   class   in  microme- 
ters,  and 

AD=dif ference   between   the   two 
classes   in  micrometers. 

The   coefficients  A.±   and   B^   are  expressed 
as 

Aj-1.136  EXP(ZS),   and  [25] 

B^-0.152   EXP(YS),  [26] 

with  Zs  and  Ys  given  by 

Zs—  6.68E-06*f-0.0903*B+1.19E-04 

*COr-3.42E-06*R0p-20400*I,   and        [27] 


Ys=3 . 28E-05*f +0 . 1 23*B-2 . 4E-04 
*COr+8.10E-06*R0p-l 1880*1, 

where  f  is  a  coefficient  and  B  is  an 
exponent  in  a  power  equation  relating 
surface  area  to  depth  expressed  as 

S  =f*Y  B 
°a  L  xd  » 


[28] 


where 

Sa=surface  area  of  the  pond  in  square 

feet,  and 
Yp=depth  in  the  pond  in  feet. 

The  coefficient  Cor  is  related  to  the 
diameter  of  the  pipe  outlet  by 

Cor=13,968*dor2,  [30] 

where 

dor=diameter  in  feet, 

ROp=volume   of    runoff   held   in   the   pond 

in  cubic   feet,   and 
I=inf iltration  rate   in  the   pond   in   feet 

per   second. 

The  equivalent  sand  diameters  for  each  of 
the  particle-size  classes  are  included  in 
table  12.   The  terms  f  and  B  are  assumed 
constant,  equal  to  7,500  and  1.5, 
respectively,  developed  from  a  range  of 
values  from  9,465  to  4,485  and  1.77  to 


[29] 
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Table  16. 

Infiltration  rates  for  each  of  the  major 

soil  texture  classes 


Soil  texture 


High  clay  soils 
Silt  soils 
High  sand  soils 
Peat  soils 


Infiltration 

rate 

(inches/h) 


0.05 
.40 
.70 

1.50 


1.10.   The  infiltration  rates  for  each  of 
the  major  soil  textures  are  shown  in 
table  16. 

The  flow  rate  leaving  the  impoundment  is 
expressed  as 


QPp=Yp°-5COr/3,600, 
where 
QpD=peak  flow  rate  leaving  the 


[31] 


<PP 


impoundment  in  cubic  feet  per  second. 


If  the  pond  is  assumed  to  be  trapezoidal 
in  shape,  the  volume  held  in  the  pond  may 
be  presented  as 

R0p=l/3(Sa*Yp).  [32] 

The  pond  depth  can  then  be  solved  as 

[33] 


Yp=(3*R0p/f)1/(B+1). 


Estimates  for  pollutant  yields  from 
feedlots  in  the  cells  are  also  calculated 
in  loop  1.   Pollutant  contributions  due 
to  these  feedlots  are  treated  as  point- 
source  pollutants  and  are  routed  into  the 
channel  with  the  others.   P  and  COD 
contributions  due  to  the  feedlots  are 
calculated  using  a  model  developed  by 
Young  et  al .  (1982b)  to  analyze  feedlot 


runoff.   The  feedlot  model  calculates  P 
and  COD  concentration  and  mass  at  both 
the  feedlot  edge  and  at  the  point  of  dis- 
charge into  a  channel  or  body  of  water. 
Concentration  at  the  feedlot  edge  is 
determined  using  the  general  relationship 


CF*R0F=C1*R01+C2*R02 , 


[34] 


where 


CF=concentration  at   the   feedlot    edge, 
R0F=volume   of    runoff   at   the    feedlot    edge, 
Ci=concentration   of    runoff   in   the 
feedlot , 
R0]^=volume   of    runoff   from  the   feedlot 
itself, 
C2=concentration  of    runoff    from  the   area 
above   the   feedlot,   and 
R02=volume   of    runoff    from  the   area   above 
the    feedlot. 

The  net   reduction  due   to    filtering   in  the 
buffer   strip   is   calculated   as 


CR=CF*(1-D1/100)*(1-D2/100), 
where 
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Cft=reduced  pollutant  concentration, 
D^=percent  reduction  of  pollutant  concen- 
tration in  overland  flow,  and 
D2=percent  reduction  of  pollutant  concen- 
tration in  any  existing  grass 
waterway. 

The  final  concentration  following  the 
buffer  is  then  calculated  as 


CT*R0T=CR*R0F+C3*R03 , 
where 
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C^final    concentration  at    the   discharge 
point , 
R0<r=total    runoff   at   the   discharge   point, 
C3=concentration  of    runoff   from   the 
area  below   the   feedlot,    and 
R03=runoff   volume    from   the   area   below   the 
feedlot. 
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Table  17. 

Background  values  for  the  total  nitrogen  concentration  in  runoff  from 

feedlots  with  unpaved  surfaces 


Location 


Total 
Nitrogen 
concentration 
(mg/L) 


Stocking 

rate 

(ft2/head) 


Reference 


Minnesota 

Ontario 

Kansas 

Nebraska 

Nebraska 

Kansas 

Ohio 

Kansas 

South  Dakota 

Kansas 

Oregon 

Wisconsin 


56 
356 
675 
854 
310 
494 

31 
364 
996 
300 
132 
486 
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212 
122 

•  •  • 

200 
200 

•  •  • 

290 
355 
570 
220 
260 
260 
500 


Young  et  al .  (1980) 
Coote  and  Hore  (1978) 
Manges  et  al .  (1975) 
Gilbertson  et  al .  (1975) 
Gilbertson  et  al.  (1970) 
Fields  (1971) 
Edwards  et  al .  (1972) 
Manges  et  al.  (1971) 
Madden  and  Dornbush  (1971) 
Miner  et  al .  (1966) 
Willrich  et  al .  (1974) 
Willrich  et  al .  (1974) 
Madison  et  al .  (1979) 


Equations  used  to  predict  the  concentra- 
tion of  total  N  in  runoff  from  feedlots 
were  developed  using  data  from  several 
feedlot  studies  done  in  the  Northern 
United  States  and  Southern  Canada. 
Background  values  for  total  N  concentra- 
tion in  runoff  from  unpaved  feedlots  were 
obtained  from  12  of  these  studies  and 
appear  in  table  17.   From  these  values,  we 
chose  300  mg/L  as  a  representative  concen- 
tration.  We  assumed  that  concentrations 
decrease  linearly  with  percent  manure  pack 
below  100  percent.   The  percent  manure 
pack  of  a  feedlot  may  be  found  by  deter- 
mining its  potential  loading.   This  can  be 
done  by  standardizing  the  potential 
pollutant  produced  by  the  animal  type  in 
the  lot  to  that  produced  by  a  1,000-pound 
beef  feeder  or  slaughter  steer. 
Comparative  factors  for  various  animal 
types  developed  using  data  from  Midwest 
Plan  Service  (1975)  and  the  American 
Society  of  Agricultural  Engineers  (1982) 


are  shown  in  table  18.   A  background 
concentration  of  12  mg/L  was  chosen  from 
the  values  in  table  19,  to  represent  the 
concentration  of  total  N  in  runoff  from 
areas  contributing  runoff  to  the  feedlot 
discharge. 

To  adjust  for  a  reduction  in  pollutant  due 
to  filtration  occurring  in  overland  flow 
once  the  runoff  leaves  the  feedlot,  the 
equation  is 


DN1=-16.8+42.3*LOG10OFT, 


[37] 


where 


DNl=percent  reduction  in  total  N  concen- 
tration due  to  overland  flow,  and 
OFT=overland  flow  time  in  seconds. 
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Table  18. 

Ratio  of  total  nitrogen  produced  by 
various  animals  to  that  produced  by  a 
1,000-pound  slaughter  steer 


Design 

weight 

Ratio 

Animal  type-*- 

(pounds) 

N 

Slaughter  steer 

1,000 

1.00 

Young  beef 

500 

.60 

Dairy  cow 

1,400 

1.68 

Young  dairy  stock 

500 

.46 

Swine 

200 

.26 

Feeder  pig 

50 

.07 

Sheep 

100 

.13 

Turkey 

10 

.02 

Chicken 

4 

.01 

Duck 

4 

.01 

Horse 

1,000 

.81 

^Data  from  Midwest  Plan  Service  (1975) 
except  for  swine,  which  is  from  American 
Society  of  Agricultural  Engineers  (1982). 


AX 


Figure  16 

Representative  watershed  cell 
with  flow  distance  along  the 
cell  noted. 


This  relationship  was  developed  using 
background  information  from  three  feedlot 
studies  (Bingham  et  al .  1978,  Dickey  and 
Vanderholm  1979,  and  Young  et  al .  1980). 
To  adjust  for  a  reduction  in  pollutant 
due  to  filtration  occurring  in  concentra- 
ted flow,  the  equation  is 


DN2=25.53+0.047*GWTC, 


[38 


where 


DN2=percent  reduction  in  total  N 

concentration  due  to  channel  flow, 
and 
GWTC=flow  time  in  the  grassed  waterway  or 
channel  in  seconds. 

This  equation  was  developed  using  data  by 
Dickey  and  Vanderholm  (1979). 


The  pollutants  from  the  feedlot  are 
treated  as  soluble  and  routed  with  the 
runoff  flow  just  as  the  other  soluble 
pollutants  are.   Once  they  reach  the 
channel,  they  are  assumed  to  remain  and 
accumulate  in  the  flow. 

Sediment  and  Nutrient  Routing 

The  sediment  routing  through  the  water- 
shed is  done  in  loops  2  and  3  of  the 
programs.   The  primary  cells  are  routed 
in  loop  2  and  the  rest  of  the  cells  in 
loop  3.   The  routing  is  done  on  a  per 
cell  and  per  particle-size  basis  proceed- 
ing from  the  headwaters  of  the  watershed 
to  its  outlet.   Figure  16  is  a  represen- 
tative cell  with  the  point  0  showing 
where  all  sediment  and  runoff  enter  the 
cell  and  point  X  showing  where  they  exit 
t  he  cell. 
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Table  19. 

Background  concentrations  for  total  nitrogen  of  runoff  waters  from 

different  land  sources 


Source 


Total 
nitrogen 
concentration 
(mg/L) 


Native  prairie — Minnesota 
Native  prairie — Oklahoma 
Continuous  corn — Minnesota 
Rotation  corn — Minnesota 
Continuous  corn — Minnesota 
Wheat — Oklahoma 
Oats — Minnesota 
Rotation  pasture — Oklahoma 
Continuous  pasture — Oklahoma 
Alf al fa — Oklahoma 
Alfalfa — Minnesota 
Fallow — Minnesota 
Aspen-birch  forest — Minnesota 


5 
9 

91 

44 

131 

2 

49 
2 
3 
4 

60 

186 

0 


Reference 


Timmons  and  Ho 
Crow  et  al.  ( 1 
Young  and  Holt 
Young  and  Holt 
Burwell  et  al . 
Olness  et  al . 
Burwell  et  al. 
Olness  et  al. 
Olness  et  al. 
Olness  et  al . 
Burwell  et  al. 
Burwell  et  al . 
Timmons  et  al . 


It  (1977) 

979) 
(1977) 
(1977) 
(1975) 

(1975) 
(1975) 

(1975) 

(1975) 

(1975) 
(1975) 
(1975) 
(1977) 


The  method  used  for  sediment  routing 
involves  equations  for  sediment  transport 
and  deposition  described  by  Foster  et  al . 
(1981)  and  Lane  (1982).   The  equation  is 
derived  from  the  steady  state  continuity 
equation. 

x 
Qs(x)=Qs(°)+QslAx/Lr-/D(x)  W  dx,    [39] 

o 
where 

Qs(x)=sediment  discharge  at  the  down- 
stream end  of  the  channel  reach  in 
pounds  per  second, 
Qs(o)=sediment  discharge  at  the  upstream 
end  of  the  channel  reach  in  pounds 
per  second, 
Qsl=lateral  sediment  inflow  rate  in 
pounds  per  second, 
x=downslope  distance  in  feet, 
Lr=reach  length  in  feet, 
D(x)=sediment  deposition  rate  at  point 

X  in  pounds  per  second-square  foot, 
and 
W=channel  width  in  feet. 


Assuming  a  constant  width,  the  deposition 
is  approximated  by 


/  D(x)Wdx=W*Ax/2*(D(o)+D(x)), 


4o; 


where 

D(o)=sediment   deposition  rate  at   the 
upstream  end   of    the   channel    in 
pounds   per  second-square   foot. 

The  deposition  is  calculated  as 

D(i)=Vss/q*(qs-gs),  [41] 

where 

D(i)=sediment  deposition  rate  at  the 

point   i   between  points  X  and   0   in 
pounds   per  second-square   foot, 


37 


VSs=particle    fall  velocity   in   feet   per 
second , 
q=runoff   rate   in   cubic   feet   per 
second-foot , 
qs=sediment    flow  rate   in  pounds    per 

second-foot,    and 
gs=effective   sediment    transport    capacity 
ia  pounds    per   second-foot. 

Effective    sediment    transport    capacity, 
gs ,    is   determined   using   a  modification   of 
the  Bagnold   (1966)   stream  power  equation: 


gs=n*k*x*vc2/Vss, 
where 
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n=an  effective  transport  factor, 
k=the  transport  capacity  factor, 
t=shear   stress   in  pounds    per   square 

foot,    and 
Vc=average   channel   velocity   in  feet   per 

second. 

The   transport    capacity   factor   is 
calculated   as 


Yi 


k=(l-eb)*es*( 


-), 
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Y  s        Yw 


where 


eb=bedload    transport   efficiency, 
es=suspended   load   transport    efficiency, 
Ys=sediraent   specific   weight   in  pounds    per 

cubic    foot,    and 
Yw=specific  weight   of   water   in  pounds   per 

cubic  foot. 


efficiency  term  was  adjusted  by  an 
effective  transport  factor,  n  (Young  et 
al.  1986).   The  value  of  n  can  be 
estimated  by 


n0.74Ef-1-98, 


[44] 


where  Ef  is  an  entrainment  function 
(Simons  and  Senturk  1976)  calculated  as 


Ef=T/((Ys-Yw)*Pd), 


[45] 


where   Pd   is    the   particle   diameter   in   feet 

The   calculated    transport    capacity 
factors,    sediment    specific  weights,    and 
particle   diameters    for   the   five 
particle-size   classes   appear   in   table   20. 
The   sediment   discharge    equation   for   each 

particle   size   is 


Qsi(x)= 


2*q(x) 


2*q(x)+AxVSsi 


Qsi(o)+Qsli 


W*Ax 


'SSI 


q(o) 


*(qSi(°)-gsi(°)) 


-Vssi/q(x)*gsi(x) 
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From  flume   studies,    the   combined 
efficiency   term,    (l-e^)*es,   has   been 
found    to   be   about   0.01    (Simons    and 
Senturk    1976).      However,    this   value   was 
determined   from  studies   using   mostly 
sand.      Since    the   actual   value    of 
(l-e^)*es   would   vary   with   the    size   of    the 
particle    being   transported,    the   combined 


where 

Qsi(x)=particle   discharge   at   the    cell 
outlet   in  pounds    per   second, 
q(x)=discharge    per  unit  width   exiting 
the   cell   in  cubic   feet   per 
second-foot , 
Ax=change    in   channel    length  across    the 
cell   in   feet, 
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Table    20. 

Sediment    specific   weights,    particle    diameters,    and    transport    capacity 

factors    for   the   5   particle-size   classes 


Specific 

Particle 

Transport 

weight. 

diameters , 

capacity 

Particle 

(lb/ft3) 

Pd   (ft) 

factor,    k 

Clay 

162.37 

6.56  *    10~6 

6.242  *   10~3 

Silt 

165.49 

3.28   *    10"5 

6.053  *    10~3 

Small   aggrega 

tes 

112.41 

1.15   *    10~4 

12.478  *    10~3 

Large    aggregates 

99.92 

1.64   *    10~3 

16.631   *    10~3 

Sand 

165.49 

6.56  *    10-4 

6.053  *    10~3 

Qs-L(o)=particle    discharge    into    the    cell    in 
pounds    per   second, 
W=average   channel  width  in   feet, 
q(o)=discharge   per  unit    of   width   into 

the   cell   in  cubic   feet   per   second- 
foot  , 

qs-^(o)=particle   discharge    per   unit    of 

width  into   the   cell   in  pounds   per 
second-foot , 

gs^(o)=particle   transport    capacity   into 

the   cell   in   pounds    per   second-foot, 

gs^(x)=particle  transport  capacity  out  of 
the  cell  in  pounds  per  second- 
foot,  and 
Qsli=lateral  particle  inflow  rate  in 
pounds  per  second. 

The  sediment  discharge  is  calculated  in 
two  periods:   the  first  period  during 
which  the  eroded  sediment  from  the  upland 
portions  of  the  cell  enter  the  channel  and 
the  remaining  period  during  which  upland 
erosion  has  stopped  but  channel  flow 
continues.   During  both  periods,  the 
sediment  flow  at  point  0  remains  constant. 
It  equals  the  total  sediment  of  that 
particular  particle  size  delivered  to  that 
point  from  all  drainage  patterns  entering 
the  cell  flowing  over  the  average  flow 
duration  calculated  at  points  0  and  X. 

The  first  period  lasts  until  overland  flow 
has  ended,  the  duration  being  equal  to  the 
time  of  overland  flow.   During  the  second 


period,  up  until  channel  flow  has  stopped, 
the  contribution  of  sediment  from  over- 
land flow  equals  0. 

The  flow  rates  in  equation  46,  q(o)  and 
q(x),  are  peak  flow  rates  divided  by  the 
appropriate  channel  width.   The  equation 
used  to  determine  the  peak  flow  rate  was 
developed  by  Smith  and  Williams  (1980)  for 
use  in  CREAMS. 

QP=8.484*A°'7*SC0.159*RF(0.824 


*A0.0166 


)(l2/(A*43560)r0-187,   [47] 


where 


QP=peak  discharge   in  cubic   feet   per 
second , 
A=drainage   area   in  acres, 
Sc=channel    slope    in   feet   per   foot, 
RF=runoff   volume    in  inches,    and 
Lc=channel    length   in   feet. 

At  point  0,  the  drainage  area  is  the 
entire  drainage  area  above  the  cell  that 
emptied  into  the  cell,  the  channel  slope 
is  that  of  the  cell  in  question,  the 
runoff  is  the  average  effective  runoff  for 
all  the  area  draining  into  the  cell,  and 
the  length  of  flow  is  the  longest  path  of 
flow  into  the  cell.   At  point  X,  the 
drainage  area  includes  the  area  of  the 
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cell  itself  plus  the  area  above,  the 
channel  slope  is  that  of  the  cell,  the 
runoff  is  the  average  effective  runoff 
including  that  of  the  cell  itself,  and  the 
length  is  the  longest  drainage  path  to  the 
bottom  of  the  cell.   For  the  primary 
cells,  the  area  is  assumed  to  be  the  cell 
area  and  the  length  is  assumed  to  be  half 
the  length  across  the  cell. 

The  effective  runoff  volume  for  each  of 
the  cells  in  the  watershed  is  also  calcu- 
lated in  loop  2.   This  volume  incorporates 
the  effects  of  impoundment  terrace  systems 
into  the  model's  runoff  calculations. 
These  terraces  will  change  the  runoff 
hydrograph  for  any  cell  with  a  terrace,  as 
well  as  that  for  any  cell  in  the  drainage 
network.   The  change  in  the  runoff  hydro- 
graph  is  incorporated  into  the  model  by 
calculating  an  effective  runoff  volume  for 
each  of  the  cells  in  the  watershed. 
Because  of  the  ponding  in  the  terrace 
system  and  the  delay  due  to  ponding,  the 
actual  peak  flow  leaving  the  cell  will 
decrease.   So  the  effective  runoff 
volume — the  volume  that  contributes  to  the 
peak  flow  leaving  the  cell — will  be  less 
for  a  cell  containing  an  impoundment 
terrace  system.   To  obtain  this  effective 
runoff  volume,  the  peak  flow  leaving  the 
terrace  system  is  added  to  the  peak  flow 
for  the  rest  of  the  cell.   The  effective 
runoff  is  then  derived  using  this  total 
peak  flow,  solving  equation  47  for  the 
runoff. 

The  length  of  the  drainage  path  is  calcu- 
lated by  summing  the  individual  drainage 
lengths  across  each  of  the  cells  that  a 
specific  drainage  path  may  cross.   This 
summation  is  made  for  each  of  the  alter 
native  paths  that  enter  the  cell  in 
question,  and  a  comparison  is  made  to 
determine  the  longest.   Flow  across  a  cell 
may  be  one  of  four  lengths,  as  shown  in 
figure  17. 

When  the  discharge  rates  are  calculated, 
duration  of  channelized  flow  is  also 


Figure    17 

Possible   drainage   paths 
across  a  cell- 
determined.      The   equation  used    to   deter- 
mine  the   channelized   flow  duration  at  both 
0   and   X   is 

D=RF*3,630*A/QP,  [48] 

where 

D=duration   in   seconds,    and 
QP=peak  discharge   in  cubic   feet   per 
second . 

To  maintain  consistent  rate  between  the  0 
and  X  points,  an  average  of  these  two 
durations  is  taken.   This  average  is  then 
used  to  determine  the  particle  discharge 
rate  into  the  cell  and  eventually  the 
mass  of  sediment  leaving  the  cell. 

The  width  of  the  channel,  in  feet,  at 
point  0  or  point  X  is  calculated  using 
the  equation 

W=2.05*z(-°-625)*(l+z2)°'125*(QP 

*n/Sc0.5)0.375>  [49] 
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where 

z=channel  sideslope  in  feet  per  foot, 

and 
n=Manning's  roughness  coefficient  for  the 

channel. 

This   equation  is   based   on   the   assumption 
that   the   channel   is   triangular  and    that 
the   flow  is   uniform.      These  widths  are 
used    to   transform  the   appropriate   terms    in 
equation  46  into   a  per  unit   of   width 
basis. 

The   lateral   sediment   flow  rate  is 
determined   by  dividing   the   amount   of 
eroded   sediment   of   each   particle-size 
class   by   the   duration  of   overland   flow. 
The  mass   of   eroded   sediment   is    found  using 
equation   1   and   the  particle-size   breakdown 
as  described    in  table    14.      The   duration  of 
overland   flow  is   defined   in  equations   4 
and   5. 

The  transport  capacity  is  calculated  as 
shown  in  equation  42.  The  shear  stress 
in  this  equation  is  solved  by  using  the 
relation 

T=Yw*Rh*Sc,  [50] 

where 

Rh=hydraulic  radius  in  feet. 

For  a  triangular  channel,  the  hydraulic 
radius  is 

Rh=y/(2*(l+z2)0.5)j  [51] 

where 

y=channel  depth  in  feet. 

By   assuming  uniform  flow,    the   velocity   in 
the   channel   can  be   calculated  using 
Manning's   equation 


Combining   these   equations    and    the    flow 
rate   equation 

Q=AC*VC,  [53] 

where 

Q=channel    flow  rate   in  cubic   feet   per 
second,   and 
Ac=cross-sectional   area  of    the   channel    in 
square   feet, 

we   can  solve   for  the   shear  stress. 

Equation  50  becomes 

t=Yw*z°*375/(2*(1+z2)0.5)0.75*Sc0.813 

*(n*QP/1.49)0*375.  [54] 

The  velocity  then  becomes 

Vc=(1.49/n)0-75*z°-25/(2*(l+z2)0*5)0-5 

*Sc0.375*qP0.25.  [55] 

The  nutrient  yield  associated  with  the 
sediment  is  calculated  using  the  total 
sediment  yield  from  each  cell.   Sediment- 
transported  nutrients  are  estimated  using 
an  equation  from  the  CREAMS  nutrient 
submodel  (Frere  et  al.  1980).   The 
equation  is 


SED-=SOIL*SED*ER*0.892, 
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Vc=(1.49/n)*Sc0.5*Rh0.667. 


[52] 


where 

SED-=N  or  P  transported  by  the  sediment 

in  pounds  per  acre, 
S0IL=N  or  P  concentration  in  the  soil, 
SED=sediment  yield  in  kilograms  per 
hectare,  and 
ER=enrichment  ratio  for  N  or  P. 

The  N  concentration  in  the  soil  is  esti- 
mated as  0.001  lb  N/lb  of  soil,  and  the  P 
as  0.0005  lb  P/lb  of  soil  (Frere  et  al . 
1980).   The  enrichment  ratio  is  calculated 
as 
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ER=a*SEDb*Tf 
where 


[57; 


a=coef f icient , 
b=exponent ,    and 
Tf=a  correction  factor  for  soil   texture. 

Factors    a  and    b   are   assumed    constant    and 
equal    to    7.4   and   -0.20,    respectively. 
The   adjustment    factors,    Tf ,    appear   in 
table    21   according   to   soil    texture. 

Supplementary  Program — DBDFL 

DBDFL  is  a  computer  program  written  to 
generate  parameter  files  for  AGNPS. 
Although  use  of  DBDFL  is  not  necessary  to 
execute  the  model,  it  greatly  eases  the 
job  of  creating  the  input  data  file  for 
the  model. 

DBDFL  was  designed  to  read  data  from  MLMIS 
and  transform  it  into  a  preliminary  data 
file  for  use  in  the  model.   DBDFL  also 
contains  options  for  entering  measured  and 
observed  parameter  values  obtained  by  the 
user.   Since  the  preliminary  data  file  is 
based  on  estimates  made  by  the  computer 
using  the  MLMIS  data,  as  many  observed 
parameters  as  possible  should  be  input 
into  the  data  file. 

Table  21. 

Adjustment  factors  (Tf)  used  to  correct 
sediment-adsorbed  nutrient  enrichment 
ratios  for  sand  and  clay  soils 


Soil  texture 


Clay 
Silt 
Sand 
Peat 


1.15 

1.00 

.85 

1.50 


Table  22. 

Parameters  assumed  constant  throughout  the 
watershed  under  default  conditions  using 
program  DBDFL 


Parameter 

Constant 

Slope  shape  factor 

1  (uniform  slope) 

Channel  sideslope 

10% 

Practice  factor 

1 

Fertilization  level 

2  (average  fertili- 

zation) 

Fertilizer  avail- 

10% 

ability  factor 

Point  source 

0 

identification 

Gully  source  level 

0 

(tons) 

Impoundment  factor 

0 

Since  not  all  of  the  variables  needed  to 
run  the  model  are  available  from  the  MLMIS 
data  file,  some  assumptions  need  to  be 
made  in  order  to  derive  the  unknown 
values .   To  create  a  preliminary  data 
file,  the  computer  makes  estimates  for  the 
unknown.   These  estimates  are  referred  to 
as  default  values.   They  simulate  charac- 
teristics for  an  average  watershed  under 
worst-case  storm  conditions.   Under  these 
conditions,  some  of  the  parameters  are 
assumed  constant  throughout  the  watershed 
(table  22). 

The  actual  parameters  used  in  the 
preliminary  data  file  are  selected  as: 

1.  Cell  number:   Cells  are  numbered 
linearly,  beginning  in  the  northwest 
corner  and  proceeding  from  west  to 
east  southward  (fig.  8). 

2.  Receiving  cell  number:   Determined  by 
the  computer  using  the  cell  aspect. 
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Table  23. 

Parameter  selection  using  MLMIS  land-use  classification 


Curve 

numbe 

r  for 

Surface 

COD 

hydro 

logic 

soil  gi 

coup 

condition 

Cropping 

Roughness 

factor 

Land  use 

A 

B 

C 

D 

constant 

factor 

coefficient 

(mg/L) 

Forested 

25 

55 

70 

77 

0.59 

0.02 

0.10 

60 

Cultivated 

77 

86 

91 

94 

.22 

.43 

.06 

120 

Water 

100 

100 

100 

100 

.01 

0 

•  •  • 

0 

Marsh 

85 

85 

85 

85 

.01 

0 

•  •  « 

25 

Urban 

77 

85 

90 

92 

.01 

.01 

.02 

110 

(Residential) 

Extractive 

30 

58 

71 

78 

.80 

.10 

.10 

80 

Pasture 

49 

69 

79 

84 

.15 

.03 

.10 

50 

Urban 

89 

92 

94 

95 

.01 

.01 

.02 

80 

( Nonres  ident  ial ) 

Transportation 

98 

98 

98 

98 

.01 

.01 

.02 

110 

5. 


6. 


7. 


Curve  number:   Determined  using  the 
land-use  (an  MLMIS  input)  and  hydro- 
logic  soil  classifications  derived 
from  MLMIS  soil  texture  information. 
Table  23  shows  the  curve  numbers  for 
each  land  use  given  in  the  MLMIS  data 
files,  and  table  24  shows  the  MLMIS 
soil  texture  groups  the  hydrologic 
soil  classifications  are  derived 
from. 

Land  slope  (percent):   Taken  directly 
from  MLMIS  data  file — wherever  a 
value  of  0  is  recorded,  the  program 
assumes  a  slope  of  0.1  percent. 


Slope  shape  factor: 
in  table  22. 


Constant  shown 


Field  slope  length:   Chosen  using  the 
field  slope  length  area  number  (fig. 
11)  and  the  land  slope  (table  25). 

Channel  slope  (percent):   Taken 
directly  from  MLMIS  data  file.   If  no 
channel  exists  in  the  cell,  the 
channel  is  represented  by  concentrated 
interrill  flow,  with  a  slope  of  one- 
half  the  land  slope.   Channel  slopes 


10. 


11, 


12, 


13. 


of    zero  are  assumed   equal    to   0.5 
percent . 

Channel   sideslope    (percent): 
Constant    shown   in   table    22. 

Manning's   roughness   coefficient: 
Determined  using   an  average   value   for 
the   land-use   category  of    the   cell 
(table   23). 

Soil  erodibility  factor:   Derived 
value  is  a  medium  for  the  range  given 
in  MLMIS  records.   Values  are  rounded 
upward  to  two  digits.   If  no  value 
exists  in  the  MLMIS  file,  one  has  to 
be  input  by  the  user. 

Cropping  factor:   Determined  using 
the  established  value  for  the  given 
land  use  (table  23). 


Practice  factor: 
table  22. 


Constant    shown   in 


Surface  condition  constant: 
Determined  using  established  value 
for  the  given  land  use  (table  23). 
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Table  24. 

Breakdown  of  MLMIS  soil  classifications  into  major  soil  textures  and 

hydrologic  soil  classes 


Soil 

Hydro 

logic  soil 

texture 

MLMIS  soil  descriptions 

classification 

Sand 

Sand 

A 

Sand  to  loamy  sand 

A 

Loamy  sand 

A 

Loamy  sand  to  sandy  loam 

B 

Sandy  loam 

B 

Sandy  loam  to  loam 

B 

Loam  to  sand  and  gravel 

B 

Silt 

Loam 

B 

Loam  to  clay  loam 

B 

Loam  to  silt  loam 

B 

Loam  to  silty  clay 

B 

Loam,  clay  loam,  sand,  and  gravel 

B 

Clay  loam 

C 

Clay  loam  to  silty  clay  loam 

C 

Silt  loam 

C 

Silty  clay  loam 

C 

Silt  loam  to  silty  clay  loam 

C 

Clay 

Silty  clay  loam  to  silty  clay 

C 

Silty  clay 

C 

Silty  clay  to  clay 
Clay 

C 
D 

Peat 


Peat 


14.  Aspect:   An  MLMIS  parameter. 

15.  Soil  texture:   MLMIS  soil  classes  are 
broken  down  into  four  soil  textures 
as  shown  in  table  24. 

16.  Fertilization  level:   Constant  shown 
in  table  22. 

17.  Fertilizer  availability  factor: 
Constant  shown  in  table  22. 

18.  Point  source  identification: 
Constant  shown  in  table  22. 


19.  Gully  source  level:   Constant  shown 
in  table  22. 

20.  Chemical  oxygen  demand  (COD)  factor: 
Determined  using  the  established 
value  for  the  given  land  use  (table 
23). 

21.  Impoundment  factor:   Constant  shown 
in  table  22. 

Using  these  methods  to  determine  the 
parameter  values,  along  with  the 


44 


Table  25. 

Slope  length  in  feet  according  to  major 

field  slope  length  area  and  land  slope 


Major  field    Slope  steepness  category  (%) 

slope  length  

area 


0-2 


3-6 


7-12 


100 
100 
200 

200 


125 
250 
250 

200 


100 
200 
200 
200 


75 
150 
150 
200 


Source:   Otterby  and  Onstad  (1978). 

information  contained  in  the  MLMIS  file, 
the  program  DBDFL  can  be  used  to  generate 
input  data  files  for  AGNPS.   Further 
discussion  of  the  use  of  DBDFL  can  be 
found  in  the  AGNPS  guide  to  model  users. 

Sensitivity  Analysis 

We  performed  a  sensitivity  analysis  on  the 
model  to  determine  the  relative  change  in 
model  output  with  respect  to  the  change  in 
inputs  and  model  variables.   After  deter- 
mining standard  input  variables  and  com- 
puting base  output  values,  we  varied  the 
input  variables  over  a  range  of  values  and 
repeated  the  computations.   These  show  how 
the  model  outputs  vary  with  changes  in 
input  values.   Results  show  how  the  model 
functions  and  how  important  each  of  the 
variables  is  in  determining  the  output. 
This  analysis  is  intended  to  aid  in 
variable  estimation. 

We  performed  the  sensitivity  analysis 
using  an  agricultural  watershed  located  in 
west-central  Minnesota.   The  watershed, 
drained  by  the  Meadow  Brook,  has  a 
drainage  area  of  about  8,680  acres  and 
outlets  into  the  Big  Stone  Lake  on  the 
Minnesota-South  Dakota  border.   The 
characteristics  of  the  watershed  for  the 
year  analyzed  are  summarized  in  table  26. 


Table  26. 

A  summary  of  the  watershed  characteristics 

and  management  practices  for  the  Meadow 

Brook  watershed  during  the  1983  growing 

season 


>13   Item 


Description 


Area 
Soil 


Land  use 


Management 
practices 


8,680  acres 

43%  loam 

38%  clay  loam 

16%  clay 

3%  silty  clay  loam 

36%  row  crops 

26%  small  grain 

3%  permanent  meadow 

6%  pasture 

2%  alfalfa 

26%  idle 

1%  farmstead 

Conventional  tillage 
Below  average  fertil- 
ization 


The  inputs  varied  in  the  sensitivity 
analysis  are  described  in  table  27.   Using 
U.S.  Department  of  Agriculture's  Soil 
Conservation  Service  records  for  the 
Meadow  Brook  watershed  during  the  1983 
growing  season,  we  established  the  base 
data  file  and  determined  the  base  output 
values.   Each  variable  was  decreased  and 
increased  by  50  and  25  percent,  while  the 
others  were  kept  constant.   In  the  case  of 
the  cell  variables,  each  of  the  cell 
values  for  the  particular  variable  being 
adjusted  was  altered  by  the  percentage 
factor.   The  output  values  described  in 
table  28  were  then  analyzed  to  determine 
their  change  in  relation  to  the  base 
values.   As  each  variable  was  varied  about 
the  base  value,  we  compared  the  mean 
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Table  27. 

A  summary  of  the  input  parameters  used  in 

the  sensitivity  analysis 


Parameter 


Rain 


EI 

CN 

LS 

FSL 

CS 

CSS 

N 

K 

C 
P 
AVAIL 


Description 


Storm  rainfall;  base  values 
for  a  25-year,  24-hour 
storm  were  assumed 

Storm  energy-intensity 

SCS   curve   number 

Land   slope 

Field  slope  length 

Channel  slope 

Channel  sideslope 

Manning's  roughness  coeffi- 
cient for  the  channel 

Soil  erodibility  factor 

Cropping  factor 

Practice  factor 

Fertilizer  availability 
factor 


Table  28. 

A  summary  of  the  output  parameters  at  the 
watershed  outlet  analyzed  in  the  sensitiv- 
ity analysis 


Parameter 


SD 


COD, 


Description 


Sediment  yield  (tons/acre) 
Nitrogen  yield  associated 

with  the  sediment 

(lb/acre) 
Soluble  nitrogen  yield 

(lb/acre) 
Phosphorus  yield  associated 

with  the  sediment 

(lb/ acre) 
Soluble  phosphorus  yield 

(lb/ acre) 
Soluble  chemical  oxygen 

demand  yield  (lb/ acre) 


output  values  with  the  means  from  the  base 
value  predictions  as  a  measure  of 
sensitivity. 

The  analysis  for  the  variables  affecting 
sediment  yield  and  sediment-associated 
nutrient  yields  is  summarized  in  table  29. 
The  first  two  columns  show  which  parame- 
ters were  varied  and  by  how  much.   The 
variables  most  significantly  affecting  the 
sediment  yield  and  the  sediment-associated 
nutrient  yields  are  the  cell  land  slope, 
the  soil  erodibility,  the  cropping  factor, 
and  the  curve  numbers.   Care  should  be 
taken  in  estimating  these,  as  well  as  the 
practice  factor,  which  also  significantly 
affects  the  output. 

We  also  analyzed  the  water  soluble 
nutrient  yield  predictions  for  sensitiv- 
ity.  This  analysis  appears  in  table  30. 
It  can  be  seen  from  table  30  that  the 


water  soluble  nutrient  yields  are  strongly 
influenced  by  the  choice  in  curve  numbers. 
Watershed  management  practices  and  crop 
covers  should  be  closely  examined  to 
derive  accurate  estimates  for  the  cell 
curve  numbers. 

From  the  sensitivity  analysis  outlined  in 
tables  29  and  30,  it  is  possible  to  deter- 
mine which  variables  need  to  be  closely 
estimated  to  make  accurate  predictions  of 
the  watershed  yields.   Although  close 
estimates  are  desirable  for  all  input 
parameters,  greater  justification  can  be 
made  for  coarsely  estimating  those  that 
least  affect  the  major  outputs  of  the 
model. 
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Table    29. 

A  summary   of    the   sensitivity   analysis    for   the    predicted   sediment 

yield   and    the   predicted   sediment-associated    nutrient    yields 


Sediment 

Nitrogen 

Phos 

pho  ru  s 

( tons/acre) 

(1W 

acre) 

(lb/ 

acre) 

Parameter^ 

Variation 

SD 

SD/base 

Ns 

Ns/base 

Ps 

Ps/base 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Base 

% 

1.00 

1.00 

3.16 

1.00 

1.58 

1.00 

Rain 

-50 

0.62 

0.62 

2.17 

0.69 

1.08 

0.69 

-25 

.85 

.85 

2.77 

.88 

1.38 

.88 

+25 

1.11 

1.11 

3.44 

1.09 

1.72 

1.09 

+50 

1.19 

1.20 

3.65 

1.15 

1.82 

1  ,15 

EI 

-50 

.50 

.50 

1.83 

.58 

.91 

.58 

-25 

.75 

.75 

2.52 

.80 

1.26 

.80 

+25 

1.25 

1.25 

3.77 

1.19 

1.89 

1.19 

+50 

1.49 

1.50 

4.36 

1.38 

2.18 

1.38 

CN 

-50 

.29 

.29 

1.16 

.37 

.58 

.37 

-25 

.70 

.70 

2.38 

.75 

1.19 

.75 

+25 

1.17 

1.17 

3.58 

1.13 

1.79 

1.13 

+50 

1.20 

1.20 

3.65 

1.16 

1.83 

1.16 

LS 

-50 

.51 

.51 

1.85 

.58 

.92 

.58 

-25 

.71 

.72 

2.42 

.76 

1.21 

.76 

+25 

1.35 

1.35 

4.02 

1.27 

2.01 

1.27 

+50 

1.68 

1.69 

4.80 

1.52 

2.40 

1.52 

FSL 

-50 

.76 

.76 

2.55 

.81 

1.27 

.81 

-25 

.89 

.89 

2.88 

.91 

1.44 

•  91 

+25 

1.09 

1.09 

3.39 

1.07 

1.69 

1.07 

+50 

1.17 

1.17 

3.58 

1.13 

1.79 

1.13 

CS 

-50 

.89 

.89 

2.88 

.91 

1.44 

.91 

-25 

.95 

.95 

3.04 

.96 

1.52 

.96 

+25 

1.03 

1.03 

3.25 

1.03 

1.62 

1.03 

+50 

1.06 

1.06 

3.32 

1.05 

1.66 

1.05 

CSS 

-50 

.76 

.77 

2.55 

.81 

1.28 

.81 

-25 

.90 

.90 

2.90 

.92 

1.45 

.92 

+25 

1.07 

1.08 

3.35 

1.06 

1.68 

1.06 

+50 

1.14 

1.14 

3.50 

1.11 

1.75 

1.11 

N 

-50 

1.16 

1.16 

3.57 

1.13 

1.78 

1.13 

-25 

1.07 

1.07 

3.33 

1.05 

1.66 

1.05 

+25 

.95 

.95 

3.04 

.96 

1.52 

.96 

+50 

.91 

.91 

2.93 

.93 

1.47 

.93 

K 

-50 

.50 

.50 

1.83 

.58 

.91 

.58 

-25 

.75 

.75 

2.52 

.80 

1.26 

.80 

+25 

1.25 

1.25 

3.77 

1.19 

1.89 

1.19 

+50 

1.49 

1.50 

4.36 

1.38 

2.18 

1.38 
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Table   29 — Continued 

A.  summary  of    the   sensitivity   analysis    for   the    predicted   sediment 

yield   and    the   predicted   sediment-associated   nutrient    yields 


Sediment 

Nitrogen 

Phos 

phorus 

Variation 

( tons/acre) 
SD     SD/base 

(lb/ 

acre) 

Ns/base 

(lb/ 

acre) 

Parameter^ 

Ps 

Ps/base 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

C 

-50 

.50 

.50 

1.83 

.58 

.91 

.58 

-25 

.75 

.75 

2.52 

.80 

1.26 

.80 

+25 

1.25 

1.25 

3.77 

1.19 

1.89 

1.19 

+50 

1.37 

1.38 

4.08 

1.29 

2.04 

1.29 

P 

-50 

.50 

.50 

1.83 

.58 

.91 

.58 

-25 

.75 

.75 

2.52 

.80 

1.26 

.80 

+25 

1.00 

1.00 

3.16 

1.00 

1.58 

1.00 

+50 

1.00 

1.00 

3.16 

1.00 

1.58 

1.00 

lEI   =   storm  energy-intensity  value,    CN  =   runoff   curve   number,   LS   =   land   slope, 
FSL   =   field  slope   length,    CS   =   channel    slope,    CSS  =   channel   sideslope,   N   = 
Manning's    roughness,    K  =   soil   erodibility    factor,    C  =   cropping   management    fac- 
tor,  and   P  =   conservation  practice   factor.      The   upper   limit   for   the    parameter   CN 
was    set   at    100;    the   upper   limit   for   the   parameters   K,    C,    and   P  was    set   at    1.00. 


Table  30. 

A  summary  of  the  sensitivity  analysis  for  the  predicted  soluble 

nutrient  yields 

Parameter*  Nitrogen  Phosphorus 

and  (lb/acre)  (lb/acre)  COD  (lb/acre) 


variation  (%)    Nq    Nq/base      Pq     p  /base      C0Dq   C0Dq/base 


Base  2.37     1.00      0.84      1.00       66.24    1.00 

Rain: 

-50 
-25 
+25 
+50 

CN: 
-50 
-25 
+25 
+50 

AVAIL: 

-50 

-25 

+25 

+50 


1.51 

.64 

.57 

.67 

2.04 

.86 

.75 

.89 

2.59 

1.09 

.90 

1.06 

2.77 

1.17 

.93 

1.10 

.07 

.03 

.02 

.03 

.36 

.15 

.11 

.13 

21.55 

9.11 

8.40 

9.96 

31.23 

13.20 

12.24 

14.50 

1.43 

.61 

.47 

.56 

1.90 

.80 

.66 

.78 

2.83 

1.20 

1.03 

1.22 

3.30 

1.39 

1.22 

1.44 

20.39 

.31 

41.99 

.63 

91.99 

1.39 

118.67 

1.79 

4.08 

.06 

26.63 

.40 

109.42 

1.65 

117.52 

1.77 

66.24 

1.00 

66.24 

1.00 

66.24 

1.00 

66.24 

1.00 

^The  upper  limit  for  the  parameters  CN  (runoff  curve  number)  and  AVAIL 
(fertilizer  availability  factor)  was  set  at  100. 
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APPENDIX  1 .   STORM  RAINFALL  VOLUME 
AND  ENERGY- INTENSITY  FREQUENCY  CURVES 
FOR  MINNESOTA 


\  2"'T 


Figure    1-1 

1-year,    24-hour    rainfall   (inches). 
Source:      U.S.    Department    of   Commerce, 
Weather   Bureau    (1961). 


Figure    1-2 

1-year,    24-hour   EI.      Source:      U.S. 
Department    of    Commerce,    Weather 
Bureau    (1961). 


Figure    1-3 

2-year,    24-hour    rainfall   (inches). 
Source:      U.S.    Department    of   Commerce, 
Weather   Bureau   (1961). 
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Figure  1-4 

2-year,  24-hour  EI.   Source:   U.S. 
Department  of  Commerce,  Weather 
Bureau  (1961). 


Figure    1-5 

5-year,    24-hour   rainfall   (inches). 
Source:      U.S.    Department    of   Commerce, 
Weather   Bureau   (1961). 


Figure    1-6 

5-year,    24-hour   EI.      Source:      U.S. 
Department   of   Commerce,    Weather 
Bureau    (1961). 
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Figure    1-7 

10-year,    24-hour   rainfall   (inches). 
Source:      U.S.    Department    of   Commerce, 
Weather  Bureau    (1961). 


Figure    1-8 

10-year,    24-hour  EI.      Source:      U.S. 
Department   of   Commerce,   Weather 
Bureau    (1961). 


Figure   1-9 

25-year,    24-hour    rainfall   (inches). 
Source:      U.S.    Department   of   Commerce, 
Weather   Bureau    (1961). 
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Figure    1-10 

15-year,    24-hour  EI.      Source:      U.S. 
Department    of   Commerce,    Weather 
Bureau   (1961). 


Figure  1-11 

50-year,  24-hour  rainfall  (inches). 
Source:   U.S.  Department  of  Commerce, 
Weather  Bureau  (1961). 


Figure  1-12 

50-year,  24-hour  EI.   Source:   U.S. 
Department  of  Commerce,  Weather 
Bureau  (1961). 
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Figure   1-13 

100-year,    24-hour   rainfall   (inches). 
Source:      U.S.    Department   of   Commerce, 
Weather  Bureau    (1961). 


Figure   1-14 

100-year,    24-hour   EI.      Source:      U.S. 
Department   of   Commerce,   Weather 
Bureau   (1961). 
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APPENDIX  2.   ESTABLISHING  A  DATA  FILE 


Background  Information 

To  demonstrate  how  to  set  up  a  data  file, 
we  have  used  an  ARS  experimental  watershed 
located  near  Treynor,  IA.   A  map  of  this 
watershed  appears  in  figure  2-1.   After 
the  map  was  obtained,  we  divided  the 
watershed  into  cells.   Because  the 
watershed  was  small,  it  was  necessary  to 
use  2.5-acre  cells  even  though  they  are 
not  normally  used  in  model  analyses. 
After  creating  the  cells,  we  numbered  them 
and  established  the  drainage  pattern  of 
each  (fig.  2-2). 

Since  the  watershed  was  field-sized,  many 
of  the  variables  were  the  same  for  all 
the  cells.   The  soil  types  of  the  water- 
shed were  predominantly  Ida  and  Manona 
silt  loams,  with  a  soil  erodibility 
factor  of  0.32,  and  hydrologic  soil  group 
B.   The  field  was  also  kept  in  continuous 
corn  farmed  on  the  contour.   The  channel 
through  the  watershed  was  a  natural 
channel  but  in  poor  condition,  with  a 
roughness  coefficient  of  0.048.   The 
particular  storm  simulated  in  this  example 
was  the  default  storm,  with  a  precipita- 
tion of  4.4  inches  and  an  energy-intensity 
value  of  56.0.   The  parameter  file  appears 
in  figure  2-3. 
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Figure  2-1 

A  contour  map  for  the  example  watershed. 
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Figure  2-2 

The  example  watershed  with  the  cells  and 

drainage  patterns. 
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Figure  2-3 

An  example  parameter  file  for  input  into 

the  model  for  the  Treynor  watershed. 
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APPENDIX  3.   MINNESOTA  LAND  MANAGEMENT 
INFORMATION  SYSTEM  (MLMIS)  DATA  FILE 
EXPLANATION 


MLMIS  Data  File  Format 


Line  1  to  end  of  file 


Column 

No. 


data  were  being  recorded.   The  other  set 
of  county  highway  maps  was  also  prepared 
by  the  Highway  Department.   This  set  of 
1936  maps  was  and  is  the  only  set  of 
statewide  maps  showing  government  lot 
numbers  from  the  PLS. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 


Row  designation 

Column  designation 

Township  number 

Range  number 

Section  number 

40-acre  parcel 

Minor  watershed  number 

Land  use/land  cover 

Forest  cover 

Soil  erodibility  code 

Land  slope  (percent) 

Soil  texture  number 

River  slope  (percent) 

360°  aspect 

River/no  river   identification 

Eight-way  aspect 


Parameter  Explanation^ 


Item  3:      Public  Land   Survey   (PLS)    Township 
Number 

Description:      This  variable   records    the 
actual    township   number   from  the   original 
Public  Land   Survey.      This  variable,    along 
with   the   other   three   PLS   variables,    is 
used   to   geographically  link  other  data 
files    to  MLMIS. 

Data   Source:      The   PLS   township   numbers 
were   recorded   by  MLMIS   staff   from  two   sets 
of    county  highway  maps.      One   set   was    the 
current   general   county   highway  maps 
prepared   by   the  Minnesota  Highway 
Department   in   the    late    1960 fs   when   the 


1    Quoted   verbatim  from   "Manual    40 — MLMIS 
Statewide   40-Acre  Data."      Land   Management 
Information  Center,   Minnesota   State 
Planning  Agency,    LL-45  Metro   Square,    7th 
and   Robert,    St.    Paul,   MN  55101. 


Coding  Procedure:   PLS  section  numbers 
were  recorded  on  coding  sheets  and 
keypunched  into  the  MLMIS  file. 

Contact  Person:   Les  Maki ,  SPA/LMIC,  (612) 
296-1208. 

Data  Classes:  PLS  township  number  is 
recorded  as  a  three-digit  number  that 
varies  from  26°  N.  to  168  N. 


Item  4:   Public  Land  Survey  (PLS)  Range 
Number 

Description:   This  variable  records  the 
actual  range  number  from  the  original 
Public  Land  Survey.   Codes  were  developed 
for  range  numbers  that  are  east  and  west 
of  the  origin  line.   This  variable,  along 
with  the  other  three  PLS  variables,  is 
used  to  geographically  link  other  data 
files  to  MLMIS. 

Data  Source:   The  PLS  range  numbers  were 
recorded  by  MLMIS  staff  from  two  sets  of 
county  highway  maps.   One  set  was  the 
current  general  county  highway  maps 
prepared  by  the  Minnesota  Highway 
Department  in  the  late  1960's  when  the 
data  were  being  recorded.   The  other  set 
of  county  highway  maps  was  also  prepared 
by  the  Highway  Department.   This  set  of 
1936  maps  was  and  is  the  only  set  of 
statewide  maps  showing  government  lot 
numbers  from  the  PLS. 

Coding  Procedure:   PLS  section  numbers 
were  recorded  on  coding  sheets  and 
keypunched  into  the  MLMIS  file. 
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Contact  Person:   Les  Maki ,  SPA/LMIC,  (612)   Item  6: 
296-1208. 


Public  Land  Survey  (PLS)  40-acre 
parcel 


Data  Classes:   PLS  range  numbers  are 
recorded  as  three-digit  numbers.   The 
first  two  digits  record  the  actual  range 
number  from  the  Public  Land  Survey.  The 
third  digit  records  east  or  west  of  the 
origin  line,  -1  for  west  and  2  for  east. 


Item  5:   Public  Land  Survey  (PLS)  Section 
Number 

Description:   This  variable  records  the 
actual  section  number  from  the  original 
Public  Land  Survey.   This  variable,  along 
with  the  other  three  PLS  variables,  is 
used  to  geographically  link  other  data 
files  to  MLMIS. 

Data  Source:   The  PLS  section  numbers  were 
recorded  by  MLMIS  staff  from  two  sets  of 
county  highway  maps.   One  set  was  the 
current  general  county  highway  maps 
prepared  by  the  Minnesota  Highway 
Department  in  the  late  1960's  when  the 
data  were  being  recorded.   The  other  set 
of  county  highway  maps  was  also  prepared 
by  the  Highway  Department.   This  set  of 
1936  maps  was  and  is  the  only  set  of 
statewide  maps  showing  government  lot 
numbers  from  the  PLS. 

Coding  Procedure:   PLS  section  numbers 
were  recorded  on  coding  sheets  and 
keypunched  into  the  MLMIS  file. 

Contact  Person:   Les  Maki,  SFM/LMIC,  (612) 
296-1208 

Data  Classes:   PLS  sections  are  recorded 
as  two-digit  numbers  ranging  from  01  to 
36. 


Description:   This  variable  records  the 
quarter-quarter  section  defined  in  the 
Public  Land  Survey.   It  records  both 
quarter-quarter  sections  from  regular 
sections  and  irregular  parcels.   This 
variable,  along  with  the  other  three  PLS 
variables,  is  used  to  geographically  link 
other  data  files  to  MLMIS. 

Data  Source:   The  PLS  section  numbers  were 
recorded  by  MLMIS  staff  from  two  sets  of 
county  highway  maps.   One  set  was  the 
current  general  county  highway  maps 
prepared  by  the  Minnesota  Highway 
Department  in  the  late  1960 's  when  the 
data  were  being  recorded.   The  other  set 
of  county  highway  maps  was  also  prepared 
by  the  Highway  Department.   This  set  of 
1936  maps  was  and  is  the  only  set  of 
statewide  maps  showing  government  lot 
numbers  from  the  PLS. 

Coding  Procedure:   PLS  section  coding 
numbers  were  recorded  on  coding  sheets  and 
keypunched  into  the  MLMIS  file. 

Contact  Person:   Les  Maki,  SPA/LMIC,  (612) 
296-2613. 

Data  Classes:   PLS  40-acre  parcels  are 
recorded  as  two-digit  numbers.   The  [16] 
codes  for  regular  parcels  are  listed 
below.   The  [16]  codes  for  irregular 
40-acre  parcels  can  be  obtained  from  the 
contact  person  listed  above. 
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Legend 
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Item  7:   Minor  Watersheds 

Contact  Person:   Earl  Nordstrand, 
SPA/LMIC,  (612)  296-1202 

Documentation  incomplete — work  in 
progress. 


Data  Source:   The  primary  information 
source  was  1:90,000  scale  air  photos  from 
black  and  white  panchromatic  film.   The 
State  south  of  St.  Cloud  was  flown  in 
April  1968  and  the  northern  half  in  April 
1969.   In  both  cases,  the  flights  were 
flown  east  to  west,  with  flight  lines 
about  10  miles  apart.  Current  county 
highway  maps  (ca.  1968),  which  at  that 
time  indicated  such  cultural  features  as 
individual  houses,  were  used  as  a 
secondary  source. 

Coding   Procedure:      Land-use /land-cover 
interpretations   were   made    at    the 
University   of   Minnesota  by   a   team   of 
students   headed   by   a   professional   photo 
interpreter.      The   data  were   coded    by   a 
recorder   as    the   photo  was    being    inter- 
preted   by  a  team  of    two.      Each   township 
was    covered   by  a  Mylar  overlay  with 
section  lines   delineated.      Each   section 
was   normally  divided    into    sixteen  40-acre 
parcels.      Actual   roads   and   field   lines 
were  used   to   identify   40-acre  parcels   on 
the   ground;    the  grid  was   used  where   no 
boundaries   occurred. 

It  is  important  to  recognize  that  the 
classification  scheme  uses  two  different 
criteria — dominant  economic  activity  (use) 
and  dominant  spatial  area  (cover) — to 
determine  the  class  of  any  40-acre  parcel. 
Spatial  dominance  was  confined  to  the 
undeveloped  (in  terms  of  structures)  land 
classes  and  generally  refers  to  cover 
type.   The  "undeveloped"  classes  are 
determined  by  dominant  economic  activity. 
Both  urban  categories  and  some  transporta- 
tion types  supercede  other  cover  types, 
and  "urban  mixed"  will  supercede  all  other 
classes  except  airports  and  railyards  in 
the  "transportation"  category. 


Item  8:   Land  Use/Land  Cover,  1969 

Description:   Nine  land-use/land-cover 
classes,  interpreted  from  high-altitude 
air  photos,  are  recorded  for  each  40-acre 
parcel  in  the  State. 


Contact  Person:   Mary  Louise  Dudding, 
SPA/LMIC,  (612)  296  2720. 
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Data  Classes ;   Data  are  coded  as  one-digit   9. 
numbers  for  the  nine  classes. 


Legend 

1 


Forested — a  forty  in  which  the 
dominant  land  use  consists  of  trees. 
To  be  considered  forested,  the  forty 
must  contain  a  scattering  of  trees 
whose  crowns  cover  at  least  10  percent 
of  the  land  area. 


2.  Cultivated — a  forty  in  which  the 
dominant  use  is  land  that  has  been 
recently  tilled  or  harvested 
mechanically. 

3.  Water — a  forty  in  which  the  dominant 
land  use  is  open  and  permanent 
water. 

4.  Marsh — a  forty  in  which  the  dominant 
land  use  consists  of  nonf ores  ted, 
shallow,  permanently  wet,  vegetated 
areas . 

5.  Urban  Residential — a  forty  containing 
five  or  more  residential  dwellings  and 
commercial  buildings. 

6.  Extractive — a  forty  in  which  the 
dominant  land  use  consists  of  the 
extraction  of  minerals,  including 
ancillary  facilities.   Examples  are 
mines,  tailing  piles,  and  gravel 
pits. 

7.  Pasture  and  Open — a  forty  on 
nonf ores  ted  land  not  used  for  any 
identifiable  purpose.   Examples  are 
grazing  land  or  abandoned  farmland. 

8.  Urban  Non-residential  or  Mixed 
Residential  Development — a  forty 
containing  at  least  one  commercial, 
industrial,  or  institutional  facility 
(including  golf  courses  and 
cemeteries)  and  possibly  containing 
residential  development. 


Transportation — a  forty  in  which  the 
dominant  land  use  consists  of 
facilities  for  the  conveyance  of 
people  or  materials. 


Item  9:   Major  Forest  Types,  1977 

Description:   The  major  forest  types  are 
groups  of  tree  species  that  commonly  grow 
in  association.   Each  variable  class 
refers  to  a  distinct  forest  type,  based  on 
standard  U.S.  Department  of  Agriculture, 
Forest  Service  (FS)  definitions  of 
"commercial"  forest  lands. 

Data  Source:   Source  data  for  this 
variable  were  collected  by  the  FS  as  part 
of  the  required  decennial  inventory  of 
Minnesota  forest  resources.   The  Resources 
Evaluation  Unit,  North  Central  Forest 
Experiment  Station,  FS,  St.  Paul,  MN, 
determined  forest  types  between  spring 
1975  and  summer  1977  by  interpreting  the 
most  recent  1:15,840-1:40,000  scale  aerial 
photographs  for  each  county.   The  currency 
of  these  photographs  ranges  between  1967 
and  1977,  depending  on  the  county.   The  FS 
recorded  forest  cover  data  on  1:24,000 
ozalid  air  photo  prints  from  the  1968-69 
high-altitude  flight.   These  prints  were 
provided  by  the  Minnesota  State  Planning 
Agency  (SPA),  and  are  permanently  on  file 
at  the  North  Central  Forest  Experiment 
Station. 

Coding  Procedure:   The  Environmental 
Planning  Division  of  SPA,  in  cooperation 
with  LMIC,  coded  forest-cover  information 
in  1977.   The  most  detailed  land  record 
maps  available,  generally  U.S.  Geological 
Survey  topographic  maps,  were  used  to 
draft  public  land  survey  section  lines  on 
to  the  air  photo  prints  described  above.  A 
40-acre  grid  was  placed  over  each  section, 
and  the  dominant  forest  type  was  recorded 
for  each  40-acre  cell  that  was  mainly 
forest  land. 
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Contact  Person:   Don  Yaeger,  SPA/ MIC,        9. 

(612)  296-2613.   For  information  on  the 

data  collection  program,  contact  Burt 

Essex,  Resources  Evaluation  Unit,  U.S. 

Department  of  Agriculture,  Forest 

Service.   North  Central  Forest  Experiment    Item  10. 

Station,    (612)   642-5282. 


Non-forest   land:      a  40-acre   parcel 
that    is   not    dominantly   covered    by 
forest    trees    or   is  within   built-up 
urban  areas. 


Surface  K-factor  for  Minnesota 
Soil  Atlas  Data 


Data  Classes:      A  one-digit    code    records 
the    six   major   forest    types    in  Minnesota. 
Unproductive    forest    land    and   nonforest 
land   have   also   been   recorded. 


Legend: 

1.  Pine:   predominantly  white  pine,  red 
pine,  or  jack  pine,  singly  or  in 
combination. 

2.  Spruce-fir:   predominantly  balsam  fir, 
white  spruce,  black  spruce,  tamarack, 
or  northern  white-cedar,  singly  or  in 
combination. 

3.  Oak:   predominantly  northern  red  oak, 
white  oak,  or  burr  oak,  singly  or  in 
combination. 

4.  Elm-ash-cot tonwood:   predominantly 
lowland  elm,  black  ash,  cot tonwood,  or 
red  maple,  singly  or  in  combination. 

5.  Maple-basswood:   predominantly  sugar 
maple,  basswood,  yellow  birch  or 
upland  elm,  and  red  maple,  singly  or 
in  combination. 

6.  Aspen-birch:   predominantly  aspen, 
balsam,  poplar,  or  paper  birch,  singly 
or  in  combination. 

7.  Unproductive  forest  land:   land 
incapable  of  producing  industrial  wood 
under  natural  conditions  because  of 
poor  site.   (Note:   Usually  poorly 
drained  swamp  conifers). 


Description:   The  soil  erodibility  factor 
(K)  in  the  Universal  Soil  Loss  Equation  is 
a  quantitative  value  experimentally 
determined  for  a  given  soil.   Generally 
speaking,  a  soil  type  becomes  less  erodi- 
ble  with  a  decrease  in  the  silt  fraction 
regardless  of  the  corresponding  increase 
in  sand  fraction  or  clay  fraction.   The 
interpretations  were  developed  for  soil 
landscape  units  (upper  5  feet  of  soil 
material  only)  within  their  respective 
georaorphic  regions. 

Primary  Data:   Surface  K-factor  is  an 
interpretation  based  on  the  consideration 
of  soil  landscape  units  and  their 
respective  geomorphic  position.   A  soil 
landscape  unit  is  a  group  of  soils 
generalized  into  a  homogeneous  unit  based 
on  subsurface  soil  texture,  surface  soil 
texture,  drainage  characteristics,  and 
surface  color.   A  geomorphic  region  is  a 
physiographic  area  defined  by  topo- 
graphic relief  and  soil  parent  material. 
For  a  more  detailed  description  of  soil 
landscape  units  and  geomorphic  regions, 
refer  to  their  Manual  40  descriptions. 

Method  of  Modification:   Professional 
interpretations  of  soils  and  geomorphic 
information  were  done  by  staff  of  the 
Department  of  Soil  Science,  University  of 
Minnesota,  in  cooperation  with  the  U.S. 
Department  of  Agriculture,  Soil 
Conservation  Service. 

Contact  Person:   Bob  Smekofski,  SPA/LMIC 
at  (612)  296-1204. 
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Item  12:   Soil  Texture  in  the  Rooting  Zone 
for  Minnesota  Soil  Atlas  Data 

Description:   Soil  texture  in  the  rooting 
zone  describes  the  dominant  presence  of 
sand,  silt,  or  clay  in  the  upper  5  feet  of 
the  soil  profile.   The  interpretations 
were  developed  for  soil  landscape  units 
within  their  respective  geomorphic 
regions. 

Primary  Data:   Soil  texture  in  the  rooting 
zone  is  an  interpretation  based  on  the 
consideration  of  soil  landscape  units  and 
their  respective  geomorphic  position.   A 
soil  landscape  unit  is  a  group  of  soils 
generalized  into  a  homogeneous  unit  based 
on  subsurface  soil  texture,  surface  soil 
texture,  drainage  characteristics,  and 
surface  color.   A  geomorphic  region  is  a 
physiographic  area  defined  by  topographic 
relief  and  soil  parent  material.   For  a 
detailed  description  of  soil  landscape 
units  and  geomorphic  regions  refer  to 
their  Manual  40  descriptions. 

Method  of  Modification:   Professional 
interpretations  of  soils  and  geomorphic 
information  were  done  by  staff  of  the 
Department  of  Soil  Science,  University  of 
Minnesota  in  cooperation  with  the  U.S. 
Department  of  Agriculture,  Soil 
Conservation  Service. 

Contact  Person:   Bob  Smekofski,  SPA/LMIC 
at  (612)  296-1204. 

Legend : 


Item  11:   Land  slope  percent 
Standard  (run/rise)  for  all  forties. 


01 

Sand 

02 

Sand  to  loamy  sand 

03 

Loam  sand 

04 

Loamy  sand  to  sandy  loam 

05 

Sandy  loam 
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06  Sandy  loam  to  loam 

0  7  Loam  to  sand  and  gravel 

08  Loam 

09  Loam  to  clay  loam 

10  Loam  to  silt  loam 

11  Loam  to  silty  clay 

12  Loam,  clay  loam,  sand  and  gravel 

13  Clay  loam 

14  Clay  loam  to  silty  clay  loam 

15  Silt  loam 

16  Silt  loam  to  silty  clay  loam 

17  Silty  clay  loam 

18  Silty  clay  loam  to  silty  clay 

19  Silty  clay 

20  Silty  clay  to  clay 

21  Peat 

22  Clay 

23  Mine  dumps 

24  Water 

Item  13:   River  Slope  (percent) 

Medium  slope  of  100-m  river  slopes  for  a 
40-acre  parcel  that  contains  rivers 

Item  14:   360°  Aspect 

Item  15:   River  or  no  River 

Nonriver  =  0,  river  =  1 

Item  16:   Aspect 

Calculated  from  vector  average  of  the 
eight  nearest  neighbors 

Legend: 


0 

Peak/dip 

1 

North 

2 

Northeast 

3 

East 

4 

Southeast 

5 

South 

6 

Southwest 

7 

West 

8 

Northwest 
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APPENDIX  4.   EXAMPLE  USE  OF  THE  SUPPLE- 
MENTARY PROGRAM — DBDFL 


Background  Information 

The  watershed  chosen  to  demonstrate  the 
program  used  to  generate  data  files  from 
MLMIS  data  is  the  same  watershed  described 
and  used  in  appendix  2.   All  options  of 
the  program  DBDFL  were  run  and  appear  in 
the  tables  of  this  appendix. 

Figure  4-1  shows  the  file  used  to  simulate 
the  MLMIS  file.   Since  this  is  not  a 
Minnesota  watershed,  no  actual  MLMIS  file 
could  be  obtained.   Thus,  many  of  the 
variables  appearing  in  the  file  are  only 
examples  to  show  what  might  be  in  a  file 
for  this  watershed.   Figure  4-2  shows  the 
edit  file  created  for  the  watershed.   The 
variables  of  this  file  and  an  explanation 
of  each  may  be  found  in  table  14.   Edited 
values  were  inserted  for  the  land-use 
code,  the  curve  number,  the  roughness 
coefficient,  the  soil  erodibility,  the 
cropping  factor,  and  the  fertilization 
level. 


Figure  4-3  shows  the  computer  interaction 
and  the  eventual  data  file  generated  when 
only  the  MLMIS  file  was  used.   In  this 
file,  as  in  the  other  files  generated,  the 
storm  characteristics  were  simulated.   To 
obtain  a  better  estimate  of  the  curve 
number,  a  hydrologic  soil  classification 
of  A  was  used. 

Figure  4-4  shows  the  results  of  using  an 
edit  file  along  with  the  MLMIS  file.   In 
figure  4-4,  all  edit  variables  were 
transferred.   The  total  edit  process 
allowed  us  to  alter  the  preliminary  data 
file  into  the  same  file  generated  manually 
in  appendix  2. 
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Figure  4-1 

An  example  MLMIS  file  for  the  demonstra- 
tion watershed.   Refer  to  "MLMIS  Data 
File  Format",  appendix  3,  for  identifica- 
tion of  columns. 
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Figure  4-2 

An  example  edit  file  for  the  demonstration 
watershed.   Refer  to  table  10  for  identi- 
fication of  columns. 


70 


TREYNOR  IOWA  DEFAULT  FILE 

2.5 

O  0 

4.4   56.0 

1 

3 

86  12.0 

1  200 

6.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

2 

5 

86   6.0 

1  250 

3.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

3 

6 

86   6.0 

1  250 

3.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

4 

7 

86   9.0 

1  200 

4.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

5 

8 

86   6.0 

1  250 

3.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

6 

9 

86   9.0 

1  200 

4.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

7 

8 

86  12.0 

1  200 

6.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

8 

12 

86   7.0 

1  200 

3.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

9 

13 

86  10.0 

1  200 

5.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

10 

11 

86   6.0 

1  250 

3.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

2 

2 

10 

0 

0 

120 

0 

11 

12 

86  12.0 

1  200 

6.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

12 

17 

86   4.0 

1  250 

2.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

13 

12 

86  10.0 

1  200 

5.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

7 

2 

2 

10 

0 

0 

120 

0 

14 

13 

86   8.0 

1  200 

4.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

7 

2 

2 

10 

0 

0 

120 

0 

15 

16 

86   5.0 

1  250 

2.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

16 

17 

86   7.0 

1  200 

3.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

17 

23 

86   5.0 

1  250 

2.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

4 

2 

2 

10 

0 

0 

120 

0 

18 

17 

86   8.0 

1  200 

4.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

7 

2 

2 

10 

0 

0 

120 

0 

19 

18 

86   6.0 

1  250 

3.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

7 

2 

2 

10 

0 

0 

120 

0 

20 

21 

86   7.0 

1  200 

3.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

21 

22 

86   6.0 

1  250 

3.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

22 

23 

86   8.0 

1  200 

4.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

23 

28 

86   7.0 

1  200 

3.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

c 

120 

0 

24 

23 

86  10.0 

1  200 

5.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

7 

2 

2 

10 

0 

0 

120 

0 

25 

26 

86  13.0 

1  150 

6.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

2 

2 

10 

0 

0 

120 

0 

26 

27 

86   6.0 

1  250 

3.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

2 

2 

10 

0 

0 

120 

0 

27 

28 

86   7.0 

1  200 

3.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

28 

33 

86   4.0 

1  250 

2.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

29 

28 

86   8.0 

1  200 

4.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

7 

2 

2 

10 

0 

0 

120 

0 

30 

31 

86  10.0 

1  200 

5.0 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

31 

32 

86  11.0 

1  200 

5.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

3 

2 

2 

10 

0 

0 

120 

0 

32 

28 

86   9.0 

1  200 

4.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

2 

2 

2 

10 

0 

0 

120 

0 

33 

34 

86   3.0 

1  250 

1.5 

5.0 

.060 

.33 

.43 

1.00 

.22 

5 

2 

2 

10 

0 

0 

120 

0 

USER  INTERACTION 

Input 

MLMIS  filename 

(ex: 

'/MLMIS)  ://DEM0 

Input 

the  name 

to  be 

assigi 

led  to  the 

output 

file 

(ex: 

(3FILEN)  :@DEM0 

Do 

you 

intend  to  use 

an  edit  fi 

le  (Y 

bs  or  Nc 

)?   NO 

Input 

the  watershed  name  (1  to 

30  characters):  TREAYNOR/ IOWA  DEFAULT 

FILE 

Input 

the  cell 

area  chosen 

(10 

or  40 

acres)  : 

2.5 

Do 

you 

wish  to 

input 

storm 

parameters  (Y 

2s  or  No) 

?  NO 

Input 

the  ma jot 

■  field 

slop 

b  len 

gth  area  i 

lumb 

er :  2 

Norma] 

close,  the  new 

data 

file 

is:  @DEM0 

Figure  4-3 

User  interaction  and  an  example  data  file 
for  program  DBDFL.   Example  showing  the 
use  of  only  MLMIS  data  file  to  generate  a 
preliminary  default  parameter  file. 
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TREYNOR  IOWA  EDITTED  FILE 

2.5 

33 

4.4   56.0 

1 

3 

75  12.0  1  200 

6.0 

5.0 

.050 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

2 

5 

75   6.0  1  250 

3.0 

5.0 

.050 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

3 

6 

75   6.0  1  250 

3.0 

5.0 

.048 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

4 

7 

75   9.0  1  200 

4.5 

5.0 

.048 

.32 

.68  1.00 

.29  5 

2  0 

10 

C 

0 

170 

0 

5 

5 

75   6.0  1  250 

3.0 

5.0 

.050 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

6 

9 

75   9.0  1  200 

4.5 

5.0 

.048 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

7 

3 

75  12.0  1  200 

6.0 

5.0 

.048 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

8 

12 

75   7.0  1  200 

3.5 

5.0 

.048 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

9 

13 

75  10.0  1  200 

5.0 

5.0 

.  048 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

10 

11 

75   6.0  1  250 

3.0 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

11 

12 

75  12.0  1  200 

6.0 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

12 

17 

75   4.0  1  250 

2.0 

5.0 

.048 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

13 

12 

75  10.0  1  200 

5.0 

5.0 

.048 

.32 

.68  1.00 

.29  7 

2  0 

10 

0 

0 

170 

0 

14 

13 

75   8.0  1  200 

4.0 

5.0 

.050 

.32 

.68  1.00 

.29  7 

2  0 

10 

0 

0 

170 

0 

15 

16 

75   5.0  1  250 

2.5 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

16 

17 

75   7.0  1  200 

3.5 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

17 

23 

75   5.0  1  250 

2.5 

5.0 

.048 

.32 

.68  1.00 

.29  4 

2  0 

10 

0 

0 

170 

0 

18 

17 

75   8.0  1  200 

4.0 

5.0 

.050 

.32 

.68  1.00 

.29  7 

2  0 

10 

0 

0 

170 

0 

19 

18 

75   6.0  1  250 

3.0 

5.0 

.050 

.32 

.68  1.00 

.29  7 

2  0 

10 

0 

0 

170 

0 

20 

21 

75   7.0  1  200 

3.5 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

21 

22 

75   6.0  1  250 

3.0 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

22 

22 

75   8.0  1  200 

4.0 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

23 

23 

75   7.0  1  200 

3.5 

5.0 

.048 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

24 

23 

75  10.0  1  200 

5.0 

5.0 

.050 

.32 

.68  1.00 

.29  7 

2  0 

10 

0 

0 

170 

0 

25 

2  6 

75  13.0  1  150 

6.5 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

26 

27 

75   6.0  1  250 

3.0 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

27 

28 

75   7.0  1  200 

3.5 

5.0 

.050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

28 

33 

75   4.0  1  250 

2.0 

5.0 

.048 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

29 

23 

75   8.0  1  200 

4.0 

5.0 

.050 

.32 

.68  1.00 

.29  7 

2  0 

10 

0 

0 

170 

0 

30 

31 

75  10.0  1  200 

5.0 

5.0 

050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

31 

32 

75  11.0  1  200 

5.5 

5.0 

050 

.32 

.68  1.00 

.29  3 

2  0 

10 

0 

0 

170 

0 

32 

23 

75   9.0  1  200 

4.5 

5.0 

050 

.32 

.68  1.00 

.29  2 

2  0 

10 

0 

0 

170 

0 

33 

34 

75   3.0  1  250 

1.5 

5.0 

048 

.32 

.68  1.00 

.29  5 

2  0 

10 

0 

0 

170 

0 

USER  INTERACTION 

Input  MLMIS  filename  (ex:  //MLMIS) 

:  //DEMO 

Input  the  name  to  be  assigned  to 

the  output  file  (ex:  @FILEN) 

: @DEMO 

Do  you  intend  to  use  ar 

i  edit  file 

(Yes 

or 

No)?   YES 

Input  the  Edit  filename 

:    (ex 

.EDITF) : 

.DEMO 

Input  the  watershed  name  (1 

to  30 

characters):  TREYN0R/I0WA  EDITTED 

FILE 

Input  the  cell  area  chosen 

[10  or 

40  acres 

):  2.5 

Do  you  wish  to  input  storm  parameters 

(Yes 

or  No)? 

NO 

Inpi 

at  the  major  field  s 

lope 

length  area  number:  2 

Normal  < 

:lose,  the  new  data  file  is:  @DEM0 

Figure  4-4 

User  interaction  and  an  example  data  file 
for  program  DBDFL.   Example  shows  the 
complete  use  of  the  edit  file  as  well  as 
the  MLMIS  data  file. 
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APPENDIX  5.   CROPPING  FACTOR  TABLES 

The   cropping   factor   tables   in   this    appendix  have   been 
reproduced    from   table    5    in  Wischmeier   and    Smith   (1978) 
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TABLE    5-A. — Approximate    soil    loss    ratios    for    cotton 


Expected    final    canopy    percent   cover:  65  80  95 

Estimated   initial   percent  cover  from  defoliation   + 

stalks    down:  30  45  60 

Practice 


Number 


Tillage   operation(s) 


Soil  loss  ratio1 


COTTON   ANNUALLY:  Percent 
] . . .  .None: 

Defoliation    to   Dec.    31  36  24  15 
Jan.    1    to    Feb.    or   Mar.   tillage: 

Cot    Rd    only  52  41  32 

Rd   &    20   percent   cover  vol    veg-  32  26  20 

Rd   &    30   percent   cover   vol    veg  26  20  14 
2  .  .      Chisel   plow    soon    offer   cot    harvest: 

Chiseling   to   Dec.   31  40  31  24 

Jan.    1    to   sprg   tillage  56  47  40 
3.  .      Fall   disk   after   chisel: 

Disking    to    Dec.    31  53  45  37 

Jan.   1   to  sprg  tillage  62  54  47 
4...    Chisel   plow   Feb~Mar,    no   prior   tillage: 

Cot  Rd   only  50  42  35 

Rd  &   20   percent   vol   veg  39  33  28 

Rd   &  30   percent  vol   veg  34  29  25 
5          Bed   ("hip")   Feb-Mar,    no   prior   tillage: 

Cot    Rd    only  100  84  70 

Rd  &  20  percent  vol   veg  78  66  56 

Rd  &   30   percent   vol   veg  68  58  50 
Split   ridges    &    plant   after    hip,    or 
Disk   &   plant  after  chisel   (SB): 

Cot    Rd    only  61  54  47 

Rd   &   20   percent  vol   veg  53  47  41 

Rd  &   30   percent  vol   veg  50  44  38 
Cropstage    1: 

Cot    Rd    only  57  50  43 

Rd   &  20  percent   vol    veg  49  43  38 

Rd  &  30  percent  vol   veg  46  41  36 

Cropstage   2  45  39  34 

Cropstage    3  40  27  17 

6.  .  .    Bed    [hip)   after    1    prior    tillage: 

Cot   Rd  only  110  96  84 

Rd    &    20    percent    veg  94  82  72 

Rd    &    30    percent    veg  90  78  68 

Split   ridges   after   hip   (SB): 

Cot    Rd   only  66  61  52 

Rd  &   20   to   30   percent   veg  61  55  49 

Cropstage    1: 

Cot    Rd    only  60  56  49 

Rd  &  20  to  30   percent   veg  56  51  46 

Cropstage   2  47  44  38 

Cropstage  3  42  30  19 

7.  .  .  .Hip   after   2   prior   tillages: 

Cot   Rd  only  116  108  98 

Rd   &   20-30  percent   veg  108  98  88 

Split    ridges    after    hip    (SB)  67  62  57 

8.  .     .Hip  after  3  or  more  tillages:  120  110  102 

Split   ridges    after    hip    (SB)  68  64  59 

9.  .  .  .  Conventiona/   moldboord   plow    and  disk: 

Fallow    period  42  39  36 

Seedbed    period  68  64  59 

Cropstage    1  63  59  55 

Cropstage   2  49  46  43 

Cropstage    3  44  32  22 
Cropstage   4   (See   practtices   1,  2,   and  3) 

COTTON   AFTER   SOD  CROP: 

For  the  first  or  second  crop  after  a  grass  or  grass-and-legume 
meadow  has  been  turnplowed,  multiply  values  given  in  the  last  five 
lines   above   by   sod   residual    factors   from   table   5-D. 

COTTON   AFTER   SOYBEANS: 

Select  values  from  above  and   multiply  by   1.25. 

See   footnotes   at   right. 


1  Alternate    procedure    for    estimating    the   soil    loss    ratios: 

The  ratios  given  above  for  cotton  are  based  on  estimates  for  re- 
ductions in  percent  cover  through  normal  winter  loss  and  by  the  succes- 
sive tillage  operations.  Research  is  underway  in  Mississippi  to  obtain 
more  accurate  residue  data  in  relation  to  tillage  practices.  This  research 
should  provide  more  accurate  soil  loss  ratios  for  cotton  within  a  few 
years. 

Where  the  reductions  in  percent  cover  by  winter  loss  and  tillage 
operations  are  small,  the  following  procedure  may  be  used  to  compute 
soil  loss  ratios  for  the  preplant  and  seedbed  periods:  Enter  figure  6  with 
the  percentage  of  the  field  surface  covered  by  residue  mulch,  move 
vertically  to  the  upper  curve,  and  read  the  mulch  factor  on  the  scale 
at  the  left.  Multiply  this  factor  by  a  factor  selected  from  the  following 
tabulation  to  credit  for  effects  of  land-use  residual,  surface  roughness 
and    porosity. 


Productivitty 
level 


No 
tillage 


Rough 
surface 


Smoothed 
surface 


High 

Medium 

Poor 


0.66 
.71 
.75 


0.50 
.54 
.58 


0.56 
.61 
.65 


Values   for  the  bedded   period  on   slopes  of   less  than    1    percent  should 
be   estimated  at   twice  the   value  computed  above  for   rough   surfaces. 

-  Rd,   crop    residue;   vol    veg,   volunteer   vegetation. 
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TABLE  5-B. — So/7  loss  ratios  for  conditions  not  evaluated 
in  table  5 

COTTON: 

See  table   5-A. 
CROPSTAGE  4  FOR  ROWCROPS: 

Stalks    broken    and    partially    standing:    Use    col.   41. 
Stalks   standing   after   hand   picking:   Col.   4L   times   1.15. 
Stalks   shredded   without  soil   tillage:    See  table  5-C. 
Fall   chisel:    Select   values  from    lines  33-62,   seedbed  column. 
CROPSTAGE  4  FOR  SMALL  GRAIN: 

See   table   5-C. 
DOUBLE   CROPPING: 

Derive   annual   C    value    by    selecting  from    table  5   the    soil    loss   per- 
centages  for   the   successive   cropstage   periods  of   each  crop. 
ESTABLISHED  MEADOW,   FULL-YEAR   PERCENTAGES: 
Grass  and  legume  mix,  3  to  5  t  hay  0.4 

Do.  2  to  3  t  hay  .6 

Do.  1    t  hay  1.0 

Sericea,  after  second   year  1.0 

Red    clover  1.5 

Alfalfa,    lespedeza,   and   second-year   sericea  2.0 

Sweetclover  2.5 

MEADOW  SEEDING  WITHOUT   NURSE  CROP: 

Determine   appropriate   lengths  of  cropstage   periods   SB,   1,  and   2  and 
apply  values  given  for  small   grain   seeding. 
PEANUTS: 

Comparison   with   soybeans    is  suggested. 
PINEAPPLES: 

Direct    data    not   available.    Tentative    values   derived    analytically    are 
available   from   the  SCS   in   Hawaii    or   the   Western   Technical   Ser- 
vice   Center   at   Portland,    Oreg.    (Reference    5). 
SORGHUM: 

Select   values   given   for  corn,   on   the  basis   of   expected   crop   residues 
and    canopy    cover. 
SUGARBEETS: 

Direct   data    not   available.    Probably    most    nearly    comparable   to    po- 
tatoes, without  the  ridging  credit. 
SUGARCANE: 

Tentative    values    available    from    sources    given    for    pineapples. 
SUMMER    FALLOW    IN    LOW-RAINFALL    AREAS,    USE    GRAIN    OR    ROW 
CROP    RESIDUES: 
The    approximate    soil    loss    percentage    after    each    successive    tillage 
operation  may  be  obtained  from  the  following  tabulation  by  esti- 
mating  the   percent    surface   cover   after   that   tillage   and   selecting 
the    column    for    the    appropriate    amount    of    initial    residue.    The 
given   values  credit   benefits   of  the   residue  mulch,   residues   mixed 
with  soil   by  tillage,  and  the   crop  system   residual. 

Percent  cover  Initial     residue    (lbs/A) 


TABLE  5-C. — Soil  loss  ratios  (percent)  for  cropstage  4 
when  stalks  are  chopped  and  distributed  without  soil 
tillage 


by  mulch 

>  4,000 

3,000 

2,000 

1,500 

90 

4 

— 

— 

— 

80 

8 

»8 

— 

— 

70 

12 

13 

'14 

— 

60 

16 

17 

'18 

'19 

50 

20 

22 

24 

>25 

40 

25 

27 

30 

32 

30 

29 

33 

37 

39 

20 

35 

39 

44 

48 

10 

47 

55 

63 

68 

Corn  or 

Sorghum 

Soyb 

sans 

Mulch 

Tilled 

Tilled 

No 

-till  in 

Grain 

cover1 

seedbed- 

No-till 

seedbed2 

corn  rd3 

Stubble1 

20 

48 

34 

60 

42 

48 

30 

37 

26 

46 

32 

37 

40 

30 

21 

38 

26 

30 

50 

22 

15 

28 

19 

22 

60 

17 

12 

21 

16 

17 

70 

12 

8 

15 

10 

12 

80 

7 

5 

9 

6 

7 

90 

4 

3 

— 

— 

4 

95 

3 

2 

— 

— 

3 

1  Part   of   a   field   surface   directly   covered   by   pieces   of   residue   mulch. 

2  This  column   applies   for   alt   systems  other  than   no-till. 

3  Cover  after  bean  harvest  may  include  an  appreciable  number  of 
stalks   carried   over  from   the   prior  corn  crop. 

4  For  grain  with  meadow  seeding,  include  meadow  growth  in  percent 
cover  and  limit  grain  period  4  to  2  mo.  Thereafter,  classify  as  estab- 
lished   meadow. 


TABLE  5-D. — Factors  to  credit  residual  effects  of  turned 
sod1 

Factor    for    cropstage    period: 

Crop  Hay  yield  

F        SB  and  1        2  3  4 

Tons 
First   year  after  mead: 

Row  crop  or  grain              3-5  0.25         0.40  0.45  0.50  0.60 

2-3  .30           .45  .50  .55  .65 

1-2  .35           .50  .55  .60  .70 
Second   year   after   mead: 

Row    crop     3.5  .70            .80  .85  .90  .95 

2-3  .75           .85  .90  .95  1.0 

1-2  .80           .90  .95  1.0  1.0 

Spring    grain    3-5  .75  .80  .85  .95 

2-3  .80  .85  .90  1.0 

1-2  .85  .90  .95  1.0 

Winter  grain    3-5  .60  .70  .85  .95 

2-3  .65  .75  .90  1.0 

1-2  —           .70  .85  .95  1.0 

1  These  factors  are  to  be  multiplied  by  the  appropriate  soil  loss  per- 
centages selected  from  table  5.  They  are  directly  applicable  for  sod- 
forming  meadows  of  at  least  1  full  year  duration,  plowed  not  more 
than   1   month  before  final   seedbed  preparation. 

When  sod  is  fall  plowed  for  spring  planting,  the  listed  values  for  all 
cropstage  periods  are  increased  by  adding  0.02  for  each  additional 
month  by  which  the  plowing  precedes  spring  seedbed  preparation.  For 
example,  September  plowing  would  precede  May  disking  by  8  months 
and  0.02(8 — 1),  or  0.14,  would  be  added  to  each  value  in  the  table.  For 
nonsod-forming  meadows,  like  sweetclover  or  lespedeza,  multiply  the 
factors  by  1.2.  When  the  computed  value  is  greater  than  1.0,  use  as  1.0. 


1  For   grain   residue   only. 

WINTER  COVER  SEEDING   IN   ROW  CROP  STUBBLE  OR   RESIDUES: 

Define   cropstage   periods   based   on   the   cover  seeding   date  and  apply 
values   from    lines    129   to    145. 
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